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ABSTRACT 



Heat transfer tests were carried out using a small enhanced tube bundle in a pool of 
R-1 13. By accurately instrumenting five tubes within the bundle, both the convective and 
nucleate boiling regions were studied in detail, with emphasis on the 'bundle' effect (i.e. 
the effect of lower tubes in operation on upper tubes within the bundle). In addition, the 
effect of surface history and pool height on nucleation site activation/deactivation was 
studied to see how this affects the overall heat transfer and in particular, the shape of the 
hysteresis loop. From the results, recommendations can be made to improve start-up 
procedures on shipboard AC systems. 
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I. INTRODUCTION 



A. BACKGROUND 

In September 1987, due to the pressing environmental issues of ozone depletion and 
global warming, an international conference was held in Montreal, Canada. This 
conference was viewed by many as a major step toward correcting the critical 
environmental problems facing the world today. Specifically, one of the problems 
addressed was chemical products damaging the protective bubble we live in. This bubble 
or thin layer is known as the ozone layer, and is what makes the earth's atmosphere 
habitable and hence unique from any other known planet Without this protective layer, the 
earth's climate would not support life as we know it today. F*rotecting this ozone layer 
from chemical destruction is paramount 

A group of chemicals that deplete the ozone layer is chlorofluorocarbons (CFCs). 
CFCs are manmade chemicals of chlorine, fluorine and carbon and were developed almost 
60 years ago by General Motors. They are unique in that they have a combination of 
desirable properties not often found in chemistry; low in toxicity , non-flammable, non- 
corrosive, non-explosive, extremely stable and compatible with many other materials. 
However, it is this extreme stability which causes problems to the ozone layer, because 
they only break down in the upper atmosphere (stratosphere) when subjected to intense 
ultraviolet radiation. This breakdown produces chlorine which has been linked to the 
depletion of the earth's protective ozone layer. A chlorine atom destroys an ozone molecule 
(O3) by reacting with it and converting it to an ordinary oxygen molecule (O 2 ) and carbon 
monoxide. More disturbing is the fact that the chlorine atom survives the ozone-destroying 
reaction and is free to repeat the process - one C!FC molecule being able to destroy 1(X),(X)0 
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ozone molecules. Unfortunately, a typical chlorine molecule can survive for up to 100 
years, so a major problem does exist. 

Armed with these facts, 24 nations (representing the United Nations Environment 
Program (UNEP)) met and signed the Montreal Protocol in 1987. They discussed 
substances that deplete the ozone layer [Ref. 1] and called for a near-term freeze on the 
production and consumption of these substances. To be more specific, it required 
production of these chemicals to be cut back to 1986 levels followed by a two-phased 
reduction culminating in cutbacks of 50% by mid- 1998; this came into effect on July 1, 
1989. In 1990, a progress meeting was held in London where UNEP delegates agreed to 
completely phase out all CFCs by the year 2000 [Ref. 2]. This was a needed, but 
worrisome change for the U.S. Navy as they use a number of different CFCs for various 
refrigeration and solvent/cleaning needs. One such need is the approximately 1850 
shipboard air conditioning (AC) plants presently using CFC-12 (in reciprocating 
compressor plants) and R-1 14 (in centrifugal compressor plants). These AC plants not 
only provide personal comfort, but cool vital electronic components,weapons systems and 
other auxiliary equipment 

To comply with the Montreal Protocol and U.S. legislation, the Mechanical Systems 
Branch/(Code 2722) at David Taylor Research Center (DTRC) is pursuing research aimed 
at eliminating shipboard use of CFCs. As mentioned by Krinsky and Noel [Ref. 3] this 
research is to be completed in three phases: 

1. In the short term, to identify suitable alternative ozone-safe chemicals to directly 
replace CFC-12 and CFC-1 14. In order to do this, the heat transfer characteristics must be 
similar to the existing refrigerants in place and hence the need for a database exists for 
current refrigerants (i.e. R-1 14, R-12, R-1 13) so that they can be compared to the new 
proposed refrigerants. 
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2. In the longer term, to research, develop and test substitute chemicals and 
alternative technologies to replace existing CFC uses. 

3. To implement new cooling system technologies into the Fleet which do not 
depend on CFCs or their replacements. One such technology under investigation at the 
Naval Postgraduate School is thermoacoustic refrigeration, which uses sound waves to 
compress and expand harmless gas molecules (such as helium) inside a tube. 

The purpose of this thesis is to supplement DTRC's research on alternatives to CFCs 
by establishing baseline nucleate pool boiling data of R-113 from a small bundle of 
enhanced tubes representing a section of a flooded evaporator . Emphasis is placed on 
surface history affecting activation/deactivation of nucleation sites, hysteresis, and the 
affect of varying evaporator pool height. The data obtained are expected to serve as 
reference data for comparison with future proposed refrigerants from bundles of both 
smooth and enhanced tube surfaces. It is felt that what is learned from this thesis about the 
characteristics of R-113 can be applied to any general refrigerant and can help the Navy 
design future AC systems using alternative chemicals. 

B . OBJECTIVES 

The objectives of this thesis are as follows: 

1. Understand in greater detail both the convection and nucleate pool boiling 
phenomena and hysteresis effects within a small enhanced tube bundle. 

2. Vary the evaporator pool height above the tube bundle to see how this effects 
nucleation site activation and general evaporator performance. 

3. Investigate the effect of surface history on nucleation site activation/deactivation. 

4. Compare data with similar studies (i.e. same fluid or tube) done at the Naval 
Postgraduate School. 
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II. LITERATURE SURVEY 



A. GENERAL INTRODUCTION 

Due to the desirability of operating many engineering devices in the nucleate boiling 
region, many types of enhanced tubes are being studied. Nishikawa et al. [Ref. 4] noted 
that there are two primary ways to increase the heat transfer performance of a boiling tube. 
The first was to decrease the wettability of the fluid in which you are immersed. They 
suggested using a Teflon coating on the boiling surface. However, this is not very 
effective for refrigerants as they are very highly wetting on any surface and fill the cavities 
of the Teflon surface. The second way proposed was to manufacture a tube with many re- 
entrant cavities, thereby trapping the vapor and keeping the nucleation sites active. These 
re-entrant cavities are designed such that if one site is activated, neighboring sites will soon 
follow and cause the entire tube to nucleate almost instantaneously due to the 
interconnecting passageways. There are now many such tubes with re-entrant cavities on 
the market 

The two main types of reboilers analyzed today are kettle reboilers and full bundle 
boilers. The difference is described in detail in Payvar [Ref. 5]. To summarize, in the 
kettle type reboiler the tube bundle only occupies the middle lower portion of the evaporator 
leaving large areas on the sides of the bundle for fluid circulation. The full bundle boilers 
occupy the whole cross section of the lower half of the evaporator shell. As can be 
reasoned from these definitions, circulation patterns are very different for the two types of 
boilers. 

The following is a brief summary of the work done on the bundle apparatus at the 
Naval Postgraduate School; 
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1 . TJ. Murphy, September 1987 [Ref. 6] 

• Designed System 

• Conducted Preliminary Operation 

2. C.L. Anderson, June 1989 [Ref. 7] 

• Tested Smooth Tubes, R-1 13 

• Tested Smooth Tubes, R-1 14/Oil 

• Tested Finned Tubes (19 fpi), R-1 14/Oil 

3. N. Akcasayar, December 1989 [Ref. 8] 

• Tested Finned Tubes (19 fpi), R-1 14/CHl 

• Tested High Flux Tubes, R-1 14/Oil 

4. H. Eraydin, December 1990 [Ref. 9] 

• Tested High Flux Tubes, R-1 14/Oil 

• Tested Turbo-B Tubes, R-1 14/Oil 

B . EXPERIMENTAL AND THEORETICAL STUDIES 

In the literature survey, no papers were found on the Turbo-B tube operating in a 
bundle or how surface history affects Turbo-B nucleation site activation/deactivation. In 
addition, no information was found on varying the pool height above the bundle to see how 
this affects nucleation site activation and general evaporator performance. These areas 
could be very important to the U.S. Navy as well as to industry and could cause the startup 
procedure of refrigeration units to be altered, thus emphasizing the need for the objectives 
of this thesis. 

In the past decade, there have been significant improvements made to tube surfaces 
and in understanding the natural convection and nucleate boiling regions for a single tube 
apparatus and several correlations have been adapted. However, the tube bundle opens a 
new arena of topics. 
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Yilmaz and Westwater [Ref. 10] studied the effect of velocity past a smooth tube and 
how it effects the heat transfer performance. They used a 6.5 mm diameter tube in R-1 13. 
They found that with an increase in velocity, the heat transfer performance of the tube was 
better. 

Cornwell and Schuller [Ref.l 1] stated that the increases in heat transfer coefficient 
found by Yilmaz and Westwater could not entirely be explained by the increase in velocity 
past the tube. Therefore, they decided to conduct a photographic study of boiling R-1 13 in 
a bundle at one atmosphere. Their study used 241 electrically heated tubes, 19.1 mm in 
diameter and arranged in a square in-line pitch of 25.4 mm, giving a pitch-to-diameter ratio 
of 1.33. They found that the bubbles leaving the lower tubes in the bundle impacted and 
caused a sliding motion around the upper tubes. Cornwell and Schuller attribute the 
increase in the upper tubes heat transfer performance on this impact and sliding motion of 
bubbles from lower tubes. Cornwell and Schuller also showed that the heat transfer 
coefficient was the highest for the top tube in the bundle and decreased as one moved down 
the bundle; the difference was as much as four times. Cornwell [Ref. 12] later found that 
in the nucleate boiling region, sliding bubbles and liquid forced convection could account 
for all the heat transfer in the top of the bundle. 

Chan and Shoukri [Ref. 13] studied the boiling characteristics of a small enhanced 
multitube bundle in a pool of R- 1 13. The rectangular bundle had a vertical pitch of 23.8 
mm and a horizontal pitch of 31.75 mm. From there studies, they concluded that the 
convective contribution is highest for tubes near the top of the bundle due to the increased 
quality and fluid velocity across these tubes. They also noted at lower heat fluxes, the 
convective effects dominate whilst for higher heat fluxes, the primary mode of heat transfCT 
was nucleate boiling. At high heat fluxes there are numerous active nucleation sites and a 
high bubble generation rate which decreased the influence of rising bubbles generated from 
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tubes below. They also found that the convective effect from a boiling tube below an upper 
tube in a bundle has less of an effect as it gets farther from the tube of interest 

A study by Bergles et al. [Ref. 14] showed that if a tube is on the verge of nucleation, 
it burst into nucleation if tapped. This test was repeated many times with the same results. 
Several other experimenters using different tubes, refrigerants, and bundle arrangements, 
have also shown the influence of tube position on heat transfer performance in the bundle. 
Wallner [Ref. 15] used a smooth tube bundle with both square and triangular geometric 
arrangements giving pitch-to-diameter ratios of 1.33 to 1.5 respectively. He found a 50 
percent increase in the heat transfer coefficient over a single tube for tubes near the top of 
the bundle. 

Fujita et al. [Ref. 16] used both smooth and enhanced tubes to study the effects of 
tube position, distribution of heat flux, and system pressure changes (0.1 to 1 MPa) within 
the bundle. They used a triangular tube arrangement with R-113. In these tests, they 
reported similar results as above with higher tubes giving improved performance. The 
enhancement in heat transfer again was believed to be caused by bubbles generated from 
the tubes below. 

Webb et al. [Ref. 17] pointed out that in air conditioning applications, the refrigerant 
typically enters the bundle with a 15 percent quality. With this entering quality into the 
bundle, it simulates the effect of lower tubes activated in the bundle and thus will improve 
the performance of the entire bundle. Furthermore, the boiling mechanism for a flooded 
evaporator is very different from the kettle type evaporator and therefore it is very hard to 
formulate and evaluate a theoretical model for a bundle. 

Leong and Cornwell [Ref. 18] looked at the same representative slice of a flooded 
evaporator that Cornwell and Schuller [Ref. 11] and Cornwell [Ref. 12] studied. They 
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concluded that the heat transfer coefficients were the highest at the top of the bundle and 
decreased in performance as you went down the bundle. 

Based on the above-mentioned studies, all experimenters seem in agreement, 
independent of the fluid or tube used in the bundle, that the presence of boiling lower tubes 
causes a convective flow in the evaporator that gives rise to an increase in the heat transfer 
coefficient on the upper tubes. 
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III. EXPERIMENTAL APPARATUS 



A. TEST APPARATUS OVERVIEW 

The whole experimental apparatus including the auxiliary equipment and the 
evaporator/condenser test apparatus used during this investigation is shown in Figure 1. 
The following is only a general description of the whole experimental apparatus. A more 
detailed look at the condenser and evaporator is provided in Section D. Further information 
about the apparatus is provided by Murphy [Ref. 6] and Anderson [Ref. 7]. 

Looking at Figure 1, the apparatus uses three closed loops. The first loop consists of 
an 8 ton refrigeration unit located outside the laboratory which is used to cool an ethylene 
glycol/water mixture. The second loop is the ethylene glycol/water mixture flowing 
through the condenser. This flow rate through the condenser was delivered by two pumps 
which could be operated independently or together; this coolant mixture condenses the 
refrigerant vapor in the condenser and maintains system pressure. Pump number one 
provided coolant flow through the four test tubes as well as to one of the auxiliary coils. 
F^imp number two provided coolant through a manifold which distributed the coolant to the 
remaining four coils in the condenser. The third loop is the evaporator/condenser itself. 
The vapor generated in the evaporator flowed upward and condensed in the condenser 
which then returned the condensate to the evaporator via gravity. 

B. AUXILIARY EQUIPMENT 

1 . 8 Ton Refrigeration Unit 

This unit was used to cool the ethylene glycol/water mixture (coolant) to the 
desired temperature needed to condense the refrigerant vapor. For the present pool boiling 
experiments, the temperature control was set to maintain the sump at 5°C. This 
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refrigeration unit had a cooling capacity of 28 kW (8 tons) with a 30 gpm pump providing 
circulation of the coolant to achieve the set temperature. To prevent damage of the unit, 
several safety features were provided: low temperature cutout, high pressure cutout, low 
flow cutout, and a hot gas bypass which allowed for continuous operation. For operation 
procedures, see Appendix E. 

2 . Ethylene GlycolAVater Mixture 

The coolant used was a 54% by weight ethylene glycol/water mixture. This 
was used to control the system pressure by circulating it through the condenser. 

3 . Pumps 

There were two pumps available to circulate the coolant from the sump through 
the condenser. Pump number one fed four test tubes and one of the secondary condenser 
coils. Pump number two was used to feed the other four secondary condenser coils. 
Pump number one provided enough cooling for tests up to five evaporator tubes. 
However, during bundle and simulation evaporator operation, both pumps were needed at 
heat fluxes above approximately 60,000 W/m^. 

4 . Flowmeters 

Four calibrated float-type flowmeters, connected to pump number one, were 
used to measure the flowrate passing through the four test condenser tubes. One additional 
flowmeter which was connected to pump number two was used to measure the total flow to 
the secondary condenser coils. In addition, each of the four secondary condenser coils had 
a globe valve to further regulate (or shut off) flow as desired. In the test tubes, the flow 
was regulated by a ball valve. 

C. DATA ACQUISITION SYSTEM/INSTRUMENTATION 

As described by Akcasayar [Ref. 8], a Hewlett Packard HP-3497A Data Acquisition 
System, HP-9125 computer and HP-701 printer were used for data acquisition, data 
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reduction and data printing respectively. Although an HP-9826 computer and HP-7470A 
plotter can be used for final graph printing, a Macintosh Classic computer using 
CricketGraph 1.3.2 was utilized. HP Basic 3.01 was used for data reduction. At the 
beginning of every set of runs, the HP-9125 computer had to be initialized by using three 
system discs if the system was turned off. As described by Anderson [Ref. 7], type-T 
copper-constantan thermocouple measurements (mvolts) were made on the HP-3497A with 
the relay multiplexer assembly equipped with thermocouple compensation. A 20 channel 
relay multiplexer card was used to measure the voltage output from voltage and amperage 
sensors. Voltage measurements were taken from separate sensors that measured the 
voltage going to the tube bundle, simulation and auxiliary heaters. The total amperage 
going through the auxiliary and simulation heaters was measured using an American 
Aerospace Control (AAC) current sensor. The currents of each instrumented tube heater 
were measured using five identical current sensors. The voltage supplied to the other active 
tubes was also measured but the current of each active tube was not (the total current for a 
pair of active tubes was measured). This was felt to be sufficient since these tubes each had 
the same power output (1000 W) as the test tube heaters and there was no apparent reason 
to monitor each active tube heat flux individually. 

Computer channel assignments for data acquisition and array assignments are given 
in Table 1 . 

D. EVAPORATOR/CONDENSER 

An overall view of the evaporator/condenser is shown in Figure 2. The evaporator 
was designed to simulate a small portion of a refrigerant flooded evaporator. Front and 
side views of the evaporator are shown in Figures 3 and 4. It was fabricated from stainless 
steel plate and formed into a short cylinder 610 mm in diameter and 241 mm long. 
Electrically heated tubes were cantilever-mounted from the back wall of the evaporator to 



11 



permit viewing along the axis of the tubes through the lower of two glass/Plexiglass 
windows mounted in the front. The glass was used on the refrigerant side in order to 
prevent crazing or cracking of the Plexiglass and the Plexiglass gave the glass added 
strength. Figure 5 is a schematic sectional view of the evaporator that shows the four kinds 
of heated tubes used in this study. The four kinds of heater tubes were as follows: 
Instrumented, Active, Auxiliary, and Simulation. Table 2 gives the power rating for the 
heaters and the number used in the evaporator. 

The electric power can be applied separately to each set of heaters using a STACO 
240V, 23.5 KVA rheostat controller shown in Figure 6. Also, the desired number of 
instrumented tubes, active tubes, auxiliary, or simulation heaters can be turned on by using 
individual circuit breakers. The four auxiliary heaters, each with a maximum rating of 4 
kW, were installed two on each side of the test bundle to maintain the liquid pool at 
saturated conditions and to provide system pressure control. The five simulation heaters, 
also capable of 4 kW each, were mounted below the test bundle in order to simulate 
additional tube rows in a larger bundle and to provide inlet vapor quality into the bottom of 
the test bundle as suggested by Webb [Ref. 17]. The auxiliary and simulation heaters had 
to be used with caution to avoid breakdown of the refrigerant 

The tube bundle consists of instrumented, active, and dummy tubes. The location of 
each tube is represented by the respective letter I, A, and D shown in Figure 5. The test 
bundle consisted of two types of heated tubes: active tubes (marked "A") which contained 
1 kW cartridge heaters, and instrumented tubes (marked"!") which, in addition to the 1 kW 
cartridge heaters, contained six wall thermocouples. In measuring boiling heat-transfer 
coefficients, great care must be exercised with the cartridge heater and temperature 
measuring instrumentation to ensure good accuracy. For example, the type of heater used 
could have an affect on heat transfer. Bergles [Ref. 19] stated that the method of heating is 
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immaterial except for the case of a.c. resistance heating which may cause bubble generation 
in phase with the power supply. However, more recently, Jung and Bergles [Ref. 20] 
who carried out extensive pool boiling data with R-1 13, concluded that the heat transfer 
coefficient of a single tube in pool boiling is not sensitive to variations in the cartridge 
heater heat flux provided that enough thermocouples are used to measure an average wall 
temperature. Prior to this investigation, the instrumented test tubes were fabricated in a 
similar way to those used by Hahne and Muller [Ref. 21] and Wanniarachchi et al. [Ref. 
22]. The exact procedure can be found in Eraydin [Ref. 9]. Figure 7 is a cross-sectional 
sketch of an instrumented tube used during this thesis, showing the construction details and 
the location of the wall thermocouples. The thermocouples were imbeded in the wall at 
different circumferential and longitudinal positions along the heated section of the tube. 
The instrumented tubes were located along the centerline of the tube bundle, forming a 
vertical in-line column. All the instrumented and active tubes were Turbo-B tubes made by 
Wolverine. These tubes were 14.15 mm in diameter and were arranged in an equilateral 
triangular pitch (i.e. centerline-to-centerline spacing) of 19.1 mm, giving a pitch-to- 
diameter ratio of 1.35. The bundle also contained a number of unheated dummy smooth 
tubes (marked "D") that were used to guide the two-phase mixture through the bundle. The 
dummy tubes were made from commercially available 15.9 mm in diameter smooth copper 
tubing. Two vertical baffle plates were made of aluminum to restrict circulation into and 
out of the bundle from the sides . A small open space approximately 5 mm in height was 
left between the bundle and the dummy tube rack. This space allowed some liquid to enter 
the bundle and replace the vapor being generated in the bundle. However, there was also a 
space below the dummy rack that provided the majority of the circulation. Thus, 
liquid/vapor circulation was vertically upward over the five instrumented test tubes with no 
net horizontal component. The dummy rack (Figure 8) used 12 tubes made of solid 
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aluminum rods each 15.9 mm in diameter and spaced 19 mm from centerline to centerline. 
This rack had a triangular pitch arrangement with vapor retaining plates on the side and was 
used for two purposes. This dummy rack was designed to collect all rising two-phase flow 
generated by the simulation heaters and direct it to the test section. The other purpose was 
to simulate the vapor passing though a larger bundle before reaching the instrumented 
tubes. The liquid-vapor mixture after passing through the bundle was separated when it 
reached the pool surface. The vapor flowed toward the condenser and the liquid was 
circulated back toward the bottom of the bundle. 

The condenser included four test tubes (each 2 m in length) in a vertical in-line 
column and five auxiliary copper coils used to regulate the pressure in the system. It was 
designed to permit independent condensation studies of small in-line tube bundles, using 
the evaporator as a source of vapor. All these tubes and coils were cooled by a refngerated 
mixture of water and ethylene glycol. Each of the test tubes and auxiliary coils could be 
operated separately using ball or globe valves respectively to regulate the desired flow rate. 
The vapor that was produced in the evaporator was guided up toward the top of the 
condenser by a vapor shroud as shown in Figure 9. The vapor was then directed down 
over the condenser tubes and coils where it condensed. Once condensed, the force of 
gravity returned the condensate back to the evaporator. Observation windows in both the 
condenser and the evaporator allowed for easy viewing of condensation and boiling. 

E. GEOMETRY OF ENHANCED SURFACE USED 

1 . Turbo-B Tube 

The Turbo-B tube, manufactured by Wolverine Tube Inc., has a enhanced 
surface geometry. The exterior boiling enhancement is made by raising integral low fins, 
cutting diagonally across these fins, and then rolling the fins to compress them to form 
mushroom-like pedestals [Ref. 23]. This process forms numerous re-entraint 
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passageways. Figure 10 shows the surface of the tube at 25 times its actual size. The tube 
is currently available in copper, cupro-nickel, and low carbon steel. The relative 
dimensions of the tube used in the present experiments are: 

Tube material - Copper 
Outside diameter = 14.15 mm 
Inside diameter = 12.7 mm 
Enhanced surface length = 203.2 mm 
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Figure 1. Schematic View of the Apparatus 
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Figure 2. Evaporator/Condenser Schematic 
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Figure 3. Front View of Evaporator 
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HIDDEN VIEW SHOWS POSITION OF RACK OF DUMMY TUBES 
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Figure 4. Side View of Evaporator 
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Figure 5. Sectional View of Evaporator Showing Tube Bundle 
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Figure 6. Photograph of 208 V, 75 A, Variable Transformers 
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Figure 7. Thermocouple Locations on an Instrumented Boiling 
Tube and Tube Section View. 
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Figure 8. Photograph of Dummy Tube Rack 
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Figure 9. Vapor flow path in Condenser 
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Figure 10. Close-up View of Turbo-B Tube Surface (25 Times) 
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TABLE 1. COMPUTER/DATA ACQUISITION ASSIGNMENT 



Thermocouple Description 


Channel 


Array in code 


Vapor 1-Top of Condenser 


00 


T(0) 


Vapor 2-Top of Condenser 


01 


T(1) 


Vapor 3-Top of Evaporator 


02 


T{2) 


Liquid 1-Top of bundle 


03 


T(3) 


Liquid 2-Top of bundle 


04 


T(4) 


Liquid 3-Bottom of bundle 


05 


T(5) 


Tube 1,No. 1 


40 


T(6) 


Tube 1,No. 2 


41 


T(7) 


Tube 1,No. 3 


42 


T{8) 


Tube 1,No. 4 


43 


T(9) 


Tube 1,No. 5 


44 


T(10) 


Tube 1,No. 6 


45 


T(11) 


Tube 2, No. 1 


46 


T(12) 


Tube 2, No. 2 


47 


T(13) 


Tube 2, No. 3 


48 


T(14) 


Tube 2, No. 4 


49 


T(15) 


Tube 2, No. 5 


50 


T(16) 


Tube 2, No. 6 


51 


T(17) 


Tube 3, No. 1 


52 


T(18) 


Tube 3, No. 2 


53 


T(19) 


Tube 3, No. 3 


54 


T(20) 


Tube 3,No. 4 


55 


T(21) 


Tube 3, No. 5 


56 


T(22) 


Tube 3,No. 6 


57 


T{23) 


Tube 4, No. 1 


58 


T{24) 


Tube 4, No. 2 


59 


T{25) 


Tube 4, No. 3 


60 


T{26) 


Tube 4, No. 4 


61 


T{27) 


Tube 4, No. 5 


62 


T{28) 


Tube 4, No. 6 


63 


T{29) 


Tube 5, No. 1 


64 


T(30) 


Tube 5, No. 2 


65 


T(31) 


Tube 5, No. 3 


66 


T(32) 


Tube 5, No. 4 


67 


T(33) 


Tube 5, No. 5 


68 


T(34) 


Tube 5, No. 6 


69 


T(35) 
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TABLE 1. COMPUTER/DATA ACQUISITION ASSIGNMENT (CONT.) 



Amperage Sensor Description 


Channel 


Array 


Tube 1 


30 


Amp(O) 


Tube 2 


31 


Amp(1 ) 


Tube 3 


32 


Amp(2) 


Tube 4 


33 


Amp(3) 


Tube 5 


34 


Amp(4) 


Active Heater Group 1 


35 


Amp(5) 


Active Heater Group 2 


36 


Amp(6) 


Active Heater Group 3 


37 


Amp(7) 


Active Heater Group 4 


38 


Amp(8) 


Active Heater Group 5 


39 


Amp(9) 


Auxiliary Heaters 


25 


Amp(1 0) 


Simulation Heaters 


26 


Amp(1 1) 



Voltage Sensor Description 


Channel 


Array 


Instrumented/Active 


27 


Volt(O) 


Simulation Heaters 


28 


Volt(1) 


Auxiliary Heaters 


29 


Volt(2) 



TABLE 2. EVAPORATOR HEATERS 



Heater Type 


Number 


Power Rating per Heater 


Instrumented Tube Heaters 


5 


1000W 


Active Tube Heaters 


1 2 


1000W 


Auxiliary Heaters 


4 


4000W 


Simulation Heaters 


5 


4000W 
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IV. EXPERIMENTAL PROCEDURES 



A . REMOVAL OF TUBE BUNDLE FROM EVAPORATOR 

Before removing the tube bundle from the evaporator, it was ensured that the 
evaporator was free of any refrigerant and was at atmospheric pressure. The tube bundle 
support block bolts were then loosened and removed. The tube bundle was then removed 
from the back of the boiler. Care was taken when removing the tube bundle as to not 
damage the thermocouple or heater wires. It was easier to disconnect the thermocouple 
wires from the data acquisition board and the heater wires from the power panel instead of 
leaving them connected to the tubes when it was necessary to work on or clean the tubes in 
the bundle. The tube bundle consisted of five instrumented heated Turbo-B tubes, 10 
active heated Turbo-B tubes, 18 dummy smooth tubes and two dummy Turbo-B tubes. 
These two dummy Turbo-B tubes were equipped with heaters and could be connected if 
desired to study various other bundle configurations (see Anderson [Ref. 7]). 

B . BUNDLE DISASSEMBLY 

When the bundle needed to be disassembled, it was ensured that there was a clean 
working space and a container to put small parts in during the procedure. The first process 
was to remove the perspex plate attached to the aluminum baffle plate by four screws. This 
plate was modified such that the eight dummy smooth tubes surrounding the instrumented 
tubes extended through the plate (with ptfe plugs which held the tubes in the proper pitch- 
to-diameter alignment). This alteration became necessary since at high heat fluxes, the 
smooth tubes (which were screwed into the tube support block) loosened from the wall and 
rattled around. The ten screws on each aluminum baffle plate were then removed. The 
plates were then pulled off the bundle. The two outside smooth dummy tubes remained on 
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the tube bundle support block as they were countersunk into the block. The three inner 
smooth tubes could be removed from the aluminum block as they were attached only by the 
screws already removed. The 12 smooth tubes were then unscrewed from the tube bundle 
support block as seen in Figure 11. Care was taken in noting the order in which all the 
tubes were removed as it saved time in reassembly. The bundle was engraved to help in 
this process. With these tubes and baffle plates removed, only the instrumented and active 
heater tubes remained. These tubes were removed by loosening the outer O-ring 
compression plate; this allowed the tubes to be pulled from the tube support block. Care 
was taken when pulling the tube from the tube support block so as to not break the 
thermocouple or heater wires. 

C. SYSTEM CLEAN-UP 

If the system was contaminated it had to be thoroughly cleaned. In order to do this, 
the entire apparatus had to be taken completely apart and cleaned in the following manner. 

After removal of the R-1 13 from the evaporator and system at atmospheric pressure, 
all electrical connections to the bundle were disconnected and the front-viewing glass 
windows were removed. The tube bundle was then be removed as described above. The 
dummy tube rack was then taken out; at all times, the system was ventilated with a fan. 

Having removed the tubes from the tube bundle, they were individually washed with 
warm water, rinsed and then wiped down with acetone. The smooth tubes were cleaned 
with Copper Brite to remove any tarnish. They were wiped down with warm water and 
acetone. The same procedure was followed for the Turbo-B tubes except they were not 
cleaned with Copper Brite. During the cleaning process, a soft bristled toothbrush was 
used to ensure the enhanced surface was cleaned properly, exercising care not to interfere 
with the tube surface. The storage tank, if thought to be contaminated also was drained, 
vented, opened and cleaned using the same procedure. 
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D. SYSTEM REFIT 



1 . Installation of tube bundle in evaporator 

Prior to installing the evaporator tubes and support block, the bundle and 
evaporator was thoroughly wiped down with acetone. Once the evaporator tubes were 
installed in the support block, and prior to tightening the tube bundle support block nuts, 
the bundle was guided into the evaporator section. After the bundle was in position, it was 
then checked that the dummy rack was still properly positioned below the bundle and that 
the vapor thermocouple positions were still 1.75 cm above the bundle. The support block 
nuts were then tightened equally on opposite sides to give equal compression of the gasket. 
If the front window was removed during maintenance, the window was put back on very 
carefully with small equal torques applied to each nut, in turn, in a clockwise direction 
around the outer ring support After the window was in place, each tube (which extended 
through the outer O-ring compression plate) was tapped forward so as to touch the front- 
viewing-window. Once completed, the O-ring compression plate was tightened. This 
compressed the O-rings (which were slightly undersized because the tubes were so close 
together) providing a seal between the tube support block and the stainless-steel backing 
plate. Also, the compression plate had grooves for the O-rings to fit into to ensure proper 
alignment and seal. 

2 . System Leakage Check 

After the system was isolated from the atmosphere and system integrity 
restored, a Seargent Welch 10 SCFM vacuum pump was turned on and the pressure in the 
system was reduced to 25 in-Hg vacuum. The system was then secured and left untouched 
for at least 10 hours to see if there was any leakage. If there was significant leakage(>l 
inHg), then the system vacuum was broken by cracking valve R-2 (Figure 2) slowly (to 
ensure minimum moisture was allowed to enter the system). This valve was left open until 
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the system was at atmospheric conditions. The system was then pressurized with air to 
approximately 5 psig through valve R-2. Large leaks were then detected by simply 
listening to the air issuing from the system. If the leaks were small, they could be found by 
applying a small amount of soap/water solution to all joints in the system. One leak that 
was found was where the sleeve had been soldered onto the Turbo-B tube. Since this was 
found to be the source of a leak, all the Turbo-B tube joints were carefully resoldered. 
After all leaks had been corrected, the system was again subjected to a vacuum for a 
minimum of 10 hours. If the increase was less than 1/2 to 1 inHg rise in a 10 hour period, 
then the system was filled with refrigerant 

E. REFRIGERANT 
1. Fill 

a. From System Storage Tank 

A refrigerant storage tank was used to store excess R-113 during 
experiments. The storage tank prevented discharge of the R-1 13 into the atmosphere and 
made experimentation less costly. To fill the evaporator with R-113 from the storage tank, 
the system pressure was equalized to the storage container pressure by opening valves R-1 
and R-7 (Figure 2). Once the pressure was equal, then valves R-6 and R-4 were opened. 
Due to the pressure difference caused by the hydrostatic head of the liquid in the storage 
container, R-113 flowed from the storage tank into the evaporator. The amount of 
refrigerant that was transferred was controlled by throttling valve R-4 to obtain the desired 
level. Once the desired amount of refrigerant had been transferred, valves R-4, R-6, R-7, 
R-8 were closed. If required, additional R-113 was later added from a storage drum at 
atmospheric pressure to the system using valve R-2 (see next section). 
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b. From Refrigerant Storage Canister 

To fill the apparatus from the refrigerant drum, the ethylene glycol/water 
mixture was first cooled to a temperature of 5°C. Both condenser pumps were put into 
operation thus allowing coolant to flow through the condenser section. This allowed for a 
differential pressure to exist between the evaporator and ambient. A clean plastic hose was 
connected to valve R-2 (Figure 2) and the other end was end put in the refrigerant drum. 
Once this differential pressure and hose hookup was completed, refrigerant transfer 
between ambient and the evaporator was possible by slowly opening valve R-2. 

2 . Removal 

a. To System Storage Tank 

For tube replacement, system maintenance or clean-up purposes, the R- 
113 was transferred to the storage tank. The ethylene glycol/water mixture temperature 
flowing through the storage tank was cooled to less than 5°C; valves R-7 and R-8 (Figure 
2) were opened and the vacuum pump was turned on to put the storage tank under vacuum. 
Once the storage tank was under a vacuum, R-113 flowed from the evaporator to the 
storage tank because of the pressure difference. This pressure difference was maintained 
by keeping the heaters on in the evaporator and the coolant flow through the storage tank. 
As the level in the evaporator came down, power to the heaters was secured to avoid 
overheating. Any remaining R-113 that could not be transferred into the storage tank was 
drained through valve R-5 and put into a drum to be used later. 

b. To Refrigerant Storage Canister 

This procedure is primarily used to change refrigerants. The system 
vacuum can be broken by slowly cracking open R-2 (Figiu'e 2) to ensure that moisture was 
not allowed to enter the system. Once the system was at atmospheric pressure, valve R-2 
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was opened wide to allow air to enter when the refrigerant was removed. The refrigerant 
could then be drained through valve R-5 and put into a drum to be used later. 

F. OPERATION 

1 . System Lightoff, Securing and Emergency Procedures 

See Appendix E 

2 . Normal Operation 

The evaporator was filled with R-1 13 to a level of approximately 10 cm above 
the top tubes in the bundle. Prior to operating the system, the 8-ton refrigeration unit was 
run for approximately an hour to reduce the ethylene-glycol/water mixture in the sump to a 
temperature of 5°C. The pressure in the evaporator was usually 14 to 16 inHg vacuum if 
the system was left secured overnight. As the sump was brought to temperature, the data 
acquisition system and computer were turned on. This allowed the temperature in the 
system to be monitored during warm-up to saturation conditions. With this and pump 
number one running, the auxiliary heaters were then used to bring the pool up to saturation 
temperature. The auxiliary heaters were set to 1 kW nominal after pool temperature had 
been reached. The system pressure was slowly raised by allowing minimal coolant to flow 
at a low rate through the condenser tubes while heating the system with the auxiliary 
heaters. Once the required saturation temperature of 47.5°C was reached, the instrumented 
tubes were switched on and set to the desired value. This was done to prevent the tubes 
from nucleating prematurely. The heat flux of the instrumented tubes was then slowly 
increased by adjusting the rheostat. For increasing heat flux, the data were taken with 
very small heat flux increments, waiting at least 5 minutes to attain steady state conditions 
at each heat flux. At all regions of the boiling curve (and especially near the onset of 
nucleate boiling), two readings were taken for each heat flux. The bundle was 
continuously monitored through the observation windows. 
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Figure 12 shows the tube bundle arrangements used during experimentation. 
Test One was with only instrumented tube number one on in the bundle. Test two was 
with instrumented tubes one and two on in the bundle. Test three was with instrumented 
tubes one, two and three on in the bundle. Test four was with instrumented tubes one 
through four on in the bundle. Test five was with all five instrumented tubes on in the 
bundle. Test six was all five instrumented tubes plus all five pairs of active heater tubes on 
in the bundle. Test seven was the same as test six but with the simulation heaters on. 

G. DATA REDUCTION PROCEDURES 

The data reduction program ’’DRP4A" was used during the experiments for 
processing the data collected. The program was written in HP Basic 3.01 and run on an 
HP-9000 series computer. The characteristics and capabilities of this software and the 
entire listing of the program is similar to that provided by Anderson [Ref. 7]. The 
following modifications were made: 

1. Correction for pool height 

2. Installation of new thermocouple at bottom of liquid pool (bundle inlet 
temperature) 

The entire program is included in Appendix F to provide continuity for future thesis 
work. 
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Figure 11. Photograph of Tube Bundle Support Block 
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Figure 12. Tube Bundle Arrangements used During Experimentation 
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Test Four Test Five Test Six Test Seven 



V. RESULTS AND DISCUSSION 



A. GENERAL COMMENTS 

The results are presented in five sections with sub-sections as appropriate. The first 
section discusses the natural convection, nucleate pool boiling phenomena and hysteresis 
effects within a small enhanced tube bundle. The second section discusses how evaporator 
pool height above the tube bundle affects the nucleation site activation and general 
evaporator performance. The third section looks at the effects of auxiliary heaters. The 
fourth section discusses the effect of surface history on nucleation site 
activation/deactivation. The final section shows comparisons of data taken during this 
thesis with previously obtained data at the Naval Postgraduate School. 

A list of data runs conducted during this investigation may be found in Appendix A. 
All data files used in this thesis use the following filename sequence. Each file is 
composed of five sets of alpha-numeric characters used to describe the experiment 

First set (1 char.) Heat flux I(Increasing heat flux) 

D(Decieasing heat flux) 

Second set (2 char.) Type of tube used TB (Turbo-B) 

Third set (2 char.) Number of Experiment 

for accounting purposes 

To give an example, the file name ITB12 means "increasing heat flux, using a Turbo- 
B tube and the 12th data run. If more detail is desired about a data run, see Appendix A. 
All plot files are similar to data files except they start with the letter "P". All tests were 
done with the auxiliary heaters at either 0, 1, or 3 kW nominal power setting (for details see 
Appendix A). 
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All the graphs show heat flux (W/m^) along the ordinate (y axis) and wall superheat 
(K) along the abscissa (x axis). The heat flux was changed through the use of a variac 
from 600 to 100,000 W/m^ for increasing or 100,000 to 600 W/m^ for decreasing runs. 
Approximately 25-30 points were taken on the increasing heat flux and about 20 points 
were taken for the decreasing heat flux. The wall superheat is the difference between the 
corrected average tube wall temperature and the liquid saturation temperature. 

To clarify the operation of the evaporator during the above tests, we used seven basic 
bundle arrangements during the test procedures as discussed and shown in Chapter IV. 

Anderson [Ref. 7] used four types of surface preparation techniques. 

1. Surface preparation A maintains the evaporator tube(s) at a heat flux of 30 kW/m^ 
for one hour, slowly reducing heat flux to a minimum value of around 1 kW/m^, and 
beginning data collection immediately at successively increasing heat fluxes. The 
motivation for this surface preparation is a widely varying dynamic load often encountered 
in an air-conditioning machine. 

2. Surface preparation B simulates an operating refrigeration system undergoing 
frequent off-cycles. A heat flux of 30 kW/m^ at saturation conditions is maintained for one 
hour. Power is then secured for 30 minutes and then data are taken at successively 
increasing heat flux steps. 

3. Surface preparation C represents the initial start-up of a refrigeration system. 
During this process, the evaporator was brought up to temperature and pressure with the 
auxiliary heaters alone to ensure the tubes did not have any active nucleation sites. Once at 
saturation pressure, the tube was turned on and the heat flux varied from 600 to 100,000 
W/m^ to develop a heat flux versus wall superheat plot This process is very interesting as 
it introduces many complexities such as boiling hysteresis. 
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4. Surface preparation D represents a continuous operating air-conditioning system. 
It is closely modelled by a continuously decreasing heat flux starting from maximum power 
of 100 kW/m^ after it has been boiled at this power for one-half hour. Considering that the 
return flow to the boiler is normally two-phase, nucleation sites are maintained active 
throughout the evaporator tube bundle. The decreasing heat flux data runs are normally the 
most readily reproducible and are used most often in measurements of tube-bundle 
performance. This surface preparation is the most widely used for reporting the 
performance of evaporator-tube surfaces. 

Three of these surface preparation techniques were used during the course of the 
present work. Surface preparation technique C and D were used exactly as stated above 
and surface preparation technique B was the same but with a modification. Surface 
preparation technique B was modified such that power to the tube heaters could be secured 
for any amount of time but the pool temperature had to be maintained at saturation 
temperature. 

1 . Repeatability of data 

Figure 13 presents similar data taken on two different days for the top tube 
operating alone in the bundle; the repeatability is seen to be very good. The arrows on the 
figure indicate the familiar hysteresis loop. The two tests show very similar temperature 
overshoot values (indicated by the horizontal line). The point of incipience (point where 
nucleate boiling first occurs) is also very similar for both tests. 

The heat-transfer coefficient, h, is often used as a measure of tube performance: 
it is found by dividing the heat flux by the wall superheat. From Figure 13, values of h in 
the nucleate boiling region range from 1540 W/m^K at low heat flux to 13,750 W/m^K at 
high heat flux. Corresponding values of h in the convection region range from 480 
W/m^K at low heat flux to 7(X) W/m^K at medium heat flux. It can be seen from Figure 
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13 that the tube performance typically improves by a factor of five once the tube starts to 
nucleate. This explains the importance of ensuring a flooded evaporator operates in the 
nucleate boiling rather than the natural convection region. 

Also shown on Figure 13 are typical uncertainty bands for the data at both high 
and low heat flux and wall superheat. It can be seen that the greatest uncertainty occurs at 
low wall superheat and heat flux; this was as expected since the values of wall superheat 
are of the order accuracy for the thermocouples. 

Figure 14 shows the limited predictions that are available for both natural 
convection and nucleate pool boiling data. These correlations were formulated for single 
smooth tubes, but represents all that is available for enhanced surfaces and tube bundles. 
This indicates the urgent need for further work in this area. For the natural convection 
region, two correlations are shown: Churchill and Chu [Ref. 24] and Churchill and Usagi 
[Ref. 25] . The experimental data for a single tube within the bundle show an appreciable 
higher heat-transfer performance than prediction. This might be expected as the tube is an 
enhanced surface. However, in the convection region, one should not expect a large 
difference in performance between different tube surfaces as they are designed primarily for 
the nucleate region. Indeed, experiments by Murphy [Ref. 6] have shown that the data for 
enhanced surfaces falls on top of smooth tube data in the convection region. 
Consequently, the discrepancy seen is more likely to be due to circulation effects within the 
bundle. 

In the nucleate boiling region, there is also a large enhancement seen in the 
experimental data when compared with the correlation of Stephan and Abdelsalam [Ref. 
26]. However, in this region, the discrepancy can be explained due to the enhanced nature 
of the surface. Circulation effects are still likely to be present, but small for a single tube 
when compared to the large number of active nucleation sites. 
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B . BUNDLE TRENDS FOR A POOL HEIGHT OF 10 CM 



Figure 15 shows the influence of tube two on tube one. One may expect the thermal 
currents from the lower tube to further enhance the upper tube. However, with this lower 
tube heated, the performance of tube one in the convection region did not change. More 
striking is the significantly poorer performance of tube two in this region for which no 
explanation can be found. It is interesting to note that tube two (which should behave as 
an isolated single tube in the convection region, i.e. feeling no effect from the tube above) 
is now in much better agreement with the natural convection correlations. 

Also seen is that the top tube nucleates at a lower heat flux. However, once both 
tubes are nucleating, the agreement is very close, with the top tube having the slightly 
better performance. Although this slight discrepancy could be due to uncertainty in the 
data, it is more likely to be due to the influence of the bubbles coming from the tube below, 
especially when one compares it with the top tube acting alone (Figure 13). 

Figure 16 shows decreasing heat flux for two tubes. Agreement is very close, with 
the top tube performing slightly better at high heat flux, but the opposite at low heat flux. 
However, as heat flux is reduced, the uncertainty in the data increases. 

With three tubes in operation (Figure 17) , the trends are much the same. One point 
of interest is the effect on lower tubes as upper tubes start to nucleate. As tube one starts 
to nucleate, tubes two and three shift significantly to the left, indicating an increase in 
heat-transfer coefficient. However, tubes two and three remain in the convection region. 
Any further movement to the left by the upper tubes causes a corresponding shift to the left 
for the tubes below. This is probably due to increased circulation effects (i.e. as the 
bubbles fi’om the upper tubes move around the evaporator, they help to agitate the thermal 
boundary layer around the tubes below). The shift to nucleation is strictly in order as one 
moves down the bundle. 
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For decreasing heat flux (Figure 18), there is a significant drop-off in performance 
for tube three, especially at a low heat flux. It may be that as one moves down the bundle, 
the greater pressure deactivates the nucleation sites at a faster rate. 

Figures 19 to 22 show four and five tubes activated for increasing and decreasing 
heat flux. Similar trends to those already reported are seen; tube one nucleates first 
followed by two, three, etc. As each successive tube nucleates, the performance of the 
remaining lower tubes is improved slightly. In the nucleate region at high heat flux, tube 
one exhibits the best performance followed by tubes two, four, five and three. This 
somewhat strange order maybe due to uncertainty or the fact that tube three is only partially 
nucleating. At the very highest heat flux, all five tubes seem to coincide. At low heat 
fluxes, tube two gives the best heat transfer; it is unclear as to why this should be. 

Because of the unexpected behavior of tube three at high heat flux, the positions of 
tubes three and five were switched to see if this behavior was a "bundle" effect or simply 
due to the particular tube itself. Figure 23 shows the results after switching for decreasing 
heat flux. New tube five (formerly tube three) now has the poorest performance, indicating 
that it is the tube rather than the bundle causing the earlier effect It may be that air pockets 
surround the thermocouple junctions of tube five. Consequently, the tubes were kept in 
this position for the remainder of the tests, with the "bad" tube in the lowest position. 

Figures 24 and 25 show the entire bundle activated for both increasing and decreasing 
heat flux while Figures 26 and 27 show the entire bundle plus simulation heaters activated 
for increasing and decreasing heat flux. These simulation heaters give the bottom of the 
bundle an inlet quality of about 20 percent, thereby simulating a real evaporator. Similar 
trends to those seen before are evident The heat flux at which incipience occurs continues 
to decrease for all tubes, indicating the extra quality moving through the bundle helps to 
trigger "early" nucleation. As this simulates part of a real evaporator, this trend is 
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encouraging, as it is more likely to be operating in the nucleate region rather than the 
convection region. For the decreasing heat flux run, the heat transfer from the lower tubes 
has also been improved due to the increased vapor passing over them. 

As already mentioned, it is difficult to make too many direct comparisons when all 
five tubes are independently instrumented. One of the tubes may have certain induced 
errors due to manufacture that another tube does not. Consequently, of perhaps more use 
is to look at what happens to the top tube only as successive tubes are activated below. 
Figure 28 does this for all seven tests (see Chapter IV: normal operation, for description of 
tests) for an increasing heat flux. In the natural convection region, the point of incipience 
occurs at increasingly lower heat fluxes as more tubes are activated. In the nucleate boiling 
region, the heat-transfer coefficient increases as more tubes are activated. Both the natural 
convection and nucleate boiling regions are well defined indicating good reliability of the 
data. For decreasing heat flux (Figure 29), the heat-transfer coefficient is again clearly 
increased as more tubes are activated below. TTiis is more evident at low heat flux than at 
high heat flux. 

C . BUNDLE PERFORMANCE WITH VARYING POOL HEIGHT 

One of the objectives of this thesis was to see the effect of pool height on heat- 
transfer performance. Pool height here is defined as the height from the free liquid surface 
down to the top of the bundle (tube one). All experiments to date had been carried out at a 
pool height of 10 cm. Two additional heights of 0 cm and 20 cm were also studied with 
the same vapor pressure within the evaporator (i.e. to see the effect of pressure head on the 
tubes in the bundle). It is understood that a true flooded evaporator has zero pool height 
(indeed the top tubes operate in the bubbly froth). 

Figures 30 to 42 show the same seven tests for a pool height of 20 cm as previously 
shown for 10 cm. The trends of all these graphs are exactly the same as for 10 cm. The 
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effect of the greater pool height is generally to delay the point of incipience. Figures 43 and 
44 show tube one only for all seven tests (i.e. to see the effect of successive activated tubes 
below) for increasing and decreasing heat flux respectively. This to shows similar trends 
to before, with the effect of successive tubes triggering earlier nucleation on the top tube as 
well as enhancing the nucleate boiling heat-transfer coefficient 

Figures 45 to 57 show the same seven tests for a 0 cm pool height and again the 
trends are very similar. However, this time, the point of incipience occurs at the lowest 
value of heat flux. Figures 58 and 59 show tube one only for all seven tests for increasing 
and decreasing heat flux respectively. The conclusions drawn from these figures are the 
same as before. 

To summarize the effect of pool height. Figures 60 and 61 show tube one (i.e. Test 
One ) at the three pool heights for increasing and decreasing heat flux. It can be seen that 
the zero pool height data nucleates first followed by 10 cm and then 20 cm. Similar trends 
were seen regardless of what test was chosen i.e. for test seven with whole bundle plus 
simulation heaters activated, all three points of incipience (for the 3 pool heights) move to a 
lower heat flux but the zero pool height data still nucleate first. Obviously, the greater 
pressure heads felt by the tubes at larger pool heights inhibit the transition to nucleate 
boiling i.e. larger wall superheats are required to overcome the increased pressure. Further 
experiments should be conducted in the future whereby the vapor pressure above the 
bundle is varied with the pool height such that the pressure that each tube "feels" at a given 
location is the same. 

D. EFFECT OF AUXILIARY HEATERS 

For all tests reported so far, the auxiliary heaters were set to IkW nominal. The 
auxiliary heaters were needed to maintain pool temperature throughout the test. They were 
only really needed at low heat fluxes or if only a few tubes in the bundle were activated. 
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For these tests, if they were not used, the evaporator could not produce enough vapor to 
keep the pressure constant even with the condenser shut down (i.e. the walls of the test 
section acted as the condenser). As a consequence, the pool temperature would drop. 
However, the necessary operation of the auxiliary heaters has an effect on the circulation 
patterns within the bundle (this has been demonstrated in the single tube facility). 
Therefore, by having these auxiliary heaters set to the same value for each test, any effect 
they may have should be a constant factor. To see the effect that the auxiliary heaters have. 
Figure 62 shows data for three tubes on during a decreasing heat flux with the auxiliary 
heaters off compared to Figure 18 with them set at 1 kW for a 10 cm pool height. 
Essentially the effect is minimal in the nucleate region at high heat fluxes. At lower heat 
fluxes, no data could be taken with the auxiliary heaters off. In the convection region 
(again at low heat flux so there is no data with the auxiliary heaters off), the effects of 
increasing the auxiliary heater power is to increase the heat-transfer coefficient (i.e. 
reducing the wall superheat). The same effect was seen for four tubes with the auxiliary 
heaters off and set to 1 kW for a 10 cm pool height; shown in Figures 63 and 20 
respectively. Figures 64 and 30 show the same trend for tube one with the auxiliary 
heaters set at 3 kW compared to 1 kW for a pool height of 20 cm. These trends are 
probably seen due to the increased circulation within the bundle at the higher auxiliary 
heater setting. Obviously, it would be preferable to do away with the auxiliary heaters 
altogether. With R-1 14, due to the lower normal boiling point, this should be possible. 

E. NUCLEATION SITE ACTIVATION/DEACTIVATION 

Nucleation site activation/deactivation is not well understood. Once a tube is 
nucleating, there are a couple of ways that the nucleation sites can die out: 

1. Reducing the heat to the tube surface thereby allowing the fluid to collapse the 
vapor bubble and re- wet the surface due to pressure effects. 
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2. Increasing the subcooling of the liquid around the tube (i.e. causing the vapor 
bubble to condense into the fluid). 

The following tests were conducted to see how quickly the nucleation sites, once 
activated, died away if the power to the tube was cut off (i.e. to see how quickly the liquid 
re- wets the surface as discussed in mechanism one above). The pool temperature was 
maintained at saturation conditions throughout the test (i.e. there was no subcooling). All 
the tests were conducted on the top tube only. The test procedure used was the modified 
surface preparation B. The power to the heater was then turned off for a known finite time 
and then reactivated. An increasing run was then carried out. This was repeated for a 
number of different times and pool heights. The auxiliary power setting was 1 kW for all 
tests. 

Figures 65 and 66 show the data for a pool height of 20 cm. The tum-off time for the 
heater varied between five seconds and one hour. It can be seen that after only 15 s, the 
nucleation sites have all died out and the tube behaves as if it had been secured overnight 
(these are the data sets listed as 14 and 15 hours). Of course, by securing overnight, 
subcooling also plays a part in deactivating the sites. 

Figure 67 shows the data for a pool height of 10 cm. The tum-off time for the heater 
varied between five seconds and one minute. It can be seen that after one minute, the 
nucleation sites have all died out and the tube behaves as if it had been secured overnight. 
Hence, it would seem that due to the reduced pressure head trying to push the liquid into 
the sites, the time needed to deactivate all the sites is significantly longer. 

Figures 68 and 69 show the data for a pool height of 0 cm. The tum-off time for the 
heater varied between 40 seconds and 8 hours (keeping pool temperature constant 
throughout). It can be seen that even for a shut-off time of 8 hours, none of the nucleation 
sites are deactivated. The pressure head is negligible here and it is mechanism two that 
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deactivates the sites. As the pool temperature was not subcooled during this particular test, 
none of the sites are therefore deactivated. Figure 70 shows a test where the apparatus was 
secured for 2.75 hours (i.e. allowed to subcool). The deactivation seen is due to the 
subcooling that has occurred. 

Similar trends for all three pool heights have been obtained for an auxiliary heater 
power setting of 3 kW. Slight differences in the times needed for total deactivation were 
noted to be 20 seconds for 10 cm vice one minute. These data are shown in Figures 71 to 
75. 



F. COMPARISONS WITH PREVIOUS NPS DATA 

Figure 76 shows a comparison between the present data for a Turbo-B bundle in R- 
113 and the data of Anderson [Ref. 7] for a smooth tube bundle in R-113. For clarity, 
only the top tube has been shown for a decreasing heat flux. It can be seen that the heat 
transfer enhancement given by the enhanced tube is three to five times for high and low 
heat fluxes respectively. 

This is also true for the overall bundle heat transfer coefficient, giving a factor of 4 
enhancement at a typical operating heat flux of 20,000 W/m^. Also shown on the figure is 
the correlation of Stephan and Abdelsalam [Ref. 26] showing a fairly good agreement with 
the smooth tube data. 

Figure 77 shows a comparison between the present data and that of Eraydin [Ref. 9] 
who used a Turbo-B bundle in R-114. It can be seen that there is very little difference 
between the two sets of data. The R-1 14 data giving the slightly better performance. This 
could be due to the slightly better thermal conductivity of R-114 at its saturation 
temperature compared to R-113 at its saturation temperature. 
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Figure 13. Performance of Tube 1 at a 10 cm Pool Height for 
Similar Data Taken on Different Days to Show Repeatability 
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Figure 14. Performance of Tube 1 at a 10 cm Pool Height Along With 
Correlations for Natural Convection and Nucleate Pool Boiling Regions 
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Figure 15. Performance Variation of Tubes 1&2 at a 10 cm Pool 
Height During a Increasing Heat Flux (ITB05) 
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Figure 16. Performance Variation of Tubes 1&2 at a 10 cm Pool 
Height During a Decreasing Heat Flux (DTB05) 
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Figure 17. Performance Variation of Tubes 1 to 3 at a 10 cm Pool 
Height During a Increasing Heat Flux (ITB06) 
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Figure 18. Performance Variation of Tubes 1 to 3 at a 10 cm Pool 
Height During a Decreasing Heat Flux (DTB06) 
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Figure 19. Performance Variation of Tubes 1 to 4 at a 10 cm Pool 
Height During a Increasing Heat Flux (ITB07) 
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Figure 20. Performance Variation of Tubes 1 to 4 at a 10 cm Pool 
Height During a Decreasing Heat Flux (DTB07) 
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Figure 21. Performance Variation of Tubes 1 to 5 at a 10 cm Pool 
Height During a Increasing Heat Flux (ITB27) 
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Figure 22. Performance Variation of Tubes 1 to 5 at a 10 cm Pool 
Height During a Decreasing Heat Flux (DTB27) 
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Figure 23. Performance Variation of Tubes 1 to 5 at a 10 cm Pool 
Height During a Decreasing Heat Flux with Tubes 3&5 Switched (DTB30) 



58 



1 00000 



10000 



C\J 

< 

E 



X 

3 



nJ 

o 

X 



1000 



100 







































































































































































































































































p 


1-1 1c 

% O 
urbo* 


















































0 

T 


B E 


3un 


dl. 












c 
































































































m ^ 


1 






































































































i»v a. 


^ 










' 




























































































■ 




A 










n 


: 
























































































OB 




















































O 


1 


1 


A 










i 


ft 








































o 




M 


■ 


§ 








o Tube 1 




























d 


? 

M 


i 














■ Tube 2 
A Tube 3 

♦ Tube 4 
X Tube 5 




















.. .... 




















































































































































































T 


r 1 

































































































































































































































.1 



1 1 0 
(Twall-Tsat)/(K) 



100 



Figure 24. Performance Variation of the Bundle at a 10 cm Pool 
Height During a Increasing Heat Flux (ITB48) 
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Figure 25. Performance Variation of the Bundle at a 10 cm Pool 
Height During a Decreasing Heat Flux (DTB48) 
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Figure 26. Performance Variation of the Bundle plus Simulation 
Heaters at a 10 cm Pool Height During a Increasing Heat Flux (ITB49) 
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Figure 27. Performance Variation of the Bundle plus Simulation 
Heaters at a 10 cm Pool Height During a Decreasing Heat Flux (DTB49) 
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Figure 28, 10 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influenced by a Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Increasing Heat Flux 
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Figure 29. 10 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influence by a Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Decreasing Heat Flux 
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Figure 30. Performance of Tube 1 at a 20 cm Pool Height 
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Figure 31. Performance Variation of Tubes 1&2 at a 20 cm Pool 
Height During a Increasing Heat Flux (ITB34) 
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Figure 32. Performance Variation of Tubes 1&2 at a 20 cm Pool 
Height During a Decreasing Heat Flux (DTB34) 



67 



1 00000 



10000 



CVJ 

< 

E 

X 

3 



nJ 

(D 

X 



1000 



100 































Or 




















































































































































































































































F 


[-112 
% Oi 
urbo- 




















































0 

T 


B E 


lun 


dl< 


B 














































































































m 


> A 




















































K 








































































































































































































































> 


















































O 






lA 
















































A 


1 


r 












































i 


<b 


i 


'Cj 

1 


r 








o Tube 1 

■ Tube 2 
A Tube 3 
































O 


r 


ft 










































■ — ■f>— 
































































































































































































L 


























































































































































































































1 


1 

































.1 



1 0 



100 



(Twall-Tsat)/{K) 



Figure 33. Performance Variation of Tubes 1 to 3 at a 20 cm Pool 
Height During a Increasing Heat Flux (ITB35) 
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Figure 34. Performance Variation of Tubes 1 to 3 at a 20 cm Pool 
Height During a Decreasing Heat Flux (DTB35) 
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Figure 35. Performance Variation of Tubes 1 to 4 at a 20 cm Pool 
Height During a Increasing Heat Flux (ITB36) 
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Figure 36. Performance Variation of Tubes 1 to 4 at a 20 cm Pool 
Height During a Decreasing Heat Flux (DTB36) 
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Figure 37. Performance Variation of Tubes 1 to 5 at a 20 cm Pool 
Height During a Increasing Heat Flux (ITB37) 
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Figure 38. Performance Variation of Tubes 1 to 5 at a 20 cm Pool 
Height During a Decreasing Heat Flux (DTB37) 
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Figure 39. Performance Variation of the Bundle at a 20 cm Pool 
Height During a Increasing Heat Flux (ITB38) 
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Figure 40, Performance Variation of the Bundle at a 20 cm Pool 
Height During a Decreasing Heat Flux (DTB38) 
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Figure 41. Performance Variation of the Bundle plus Simulation 
Heaters at a 20 cm Pool Height During a Increasing Heat Flux (ITB40) 
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Figure 42. Performance Variation of the Bundle plus Simulation 
Heaters at a 20 cm Pool Height During a Decreasing Heat Flux (DTB40) 
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Pc'^rormance Variation of Tube 1 in a Pool of R- 
113 When InHuenced by an Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Increasing Heat Flux 
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Figure 44. 20 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influenced by an Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Decreasing Heat Flux. 
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Figure 45. Performance of Tube 1 at a 0 cm Pool Height 
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Figure 46. Performance Variation of Tubes 1&2 at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB52) 
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Figure 47. Performance Variation of Tubes 1&2 at a 0 cm Pool 
Height During a Decreasing Heat Flux (DTB52) 
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Figure 48. Performance Variation of Tubes 1 to 3 at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB53) 
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Figure 49, Performance Variation of Tubes 1 to 3 at a 0 cm Pool 
Height During a Decreasing Heat Flux (DTB53) 
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Figure 50. Performance Variation of Tubes 1 to 4 at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB55) 
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Figure 51* Performance Variation of Tubes 1 to 4 at a 0 cm Pool 
Height During a Decreasing Heat Flux (DTB55) 
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Figure 52. Performance Variation of Tubes 1 to 5 at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB56) 
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Figure 53. Performance Variation of Tubes 1 to 5 at a 0 cm Pool 
Height During a Decreasing Heat Flux (DTB56) 
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Figure 54. Performance Variation of the Bundle at a 0 cm Pool 
Height During a Increasing Heat Flux (ITB57) 



89 



100000 



10000 



C\J 

< 

E 



X 

3 



(T5 

(D 

X 



1000 



100 





























ET 










































































































































































































































R 


-112 
% O 
urbo- 


















































0 

T 


B E 


}un 


dll 












X 












































..Mi 




















































n 









































































































































































































































































































































































i 


► * 


















































a 


L 






















o Tube 1 

■ Tube 2 
A Tube 3 

♦ Tube 4 
X Tube 5 
















* 


■ 


o 

: 


4 













































































































































































































































































































































































































































































.1 



1 0 



100 



(Twall-Tsat)/(K) 



Figure 55. Performance Variation of the Bundle at a 0 cm Pool 
Height During a Decreasing Heat Flux (DTB57) 
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Figure 56. Performance Variation of the Bundle plus Simulation 
Heaters at a 0 cm Pool Height During a Increasing Heat Flux (ITB58) 
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Figure 57. Performance Variation of the Bundle plus Simulation 
Heaters at a 0 cm Pool Height During a Decreasing Heat Flux (DTB58) 



92 



100000 



10000 



CVJ 

< 

E 

§ 

X 

13 



0} 

CD 

X 



1000 



100 



.1 



R-1 13*^ 

•0% Oil- 
Turbo-B Bundle 






dUi 



ID 



it 









r 






1 10 
(Twall-Tsat)/(K) 



□ 1T(ITB51) 



♦ 2T(ITB52) 



3T(ITB53) 



♦ 4T(ITB55) 



5T(ITB56) 



o B(ITB57) 



* B+S(ITB58) 



100 



Figure 58. 0 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influenced by an Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Increasing Heat Flux 



93 



100000 



10000 



CVJ 

< 

E 

— 

5 

X 

3 



C0 

a> 

X 



1000 



100 































































































































Q 














































_ _ no/ r\il 






























Ti 


/O V^ll 

jrbo- 


B E 


un 


die 










































n 


* P 
































r 


F " 

l_ 






































































































































































































































M 






Q 


























rv-i 






























i 


R 




1: 


























■ ■ ■ ■ ■ ^ 


1 






































































































































li 


J 





































































































































































Q 1T(DTB51) 



♦ 2T(DTB52) 



3T(DTB53) 



♦ 4T(DTB55) 



5T(DTB56) 



D B(DTB57) 



A B+S(DTB58) 



1 10 
(Twall-Tsat)/(K) 



100 



Figure 59. 0 cm. Performance Variation of Tube 1 in a Pool of R- 
113 When Influenced by an Increasing Number of Heated Tubes in a 
Turbo-B Bundle During a Decreasing Heat Flux 
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Figure 60. Comparison of Tube 1 in a Pool of R-113 for Three 
Different Pool Heights During a Increasing Heat Flux 
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Figure 61. Comparison of Tube 1 in a Pool of R-113 for Three 
Different Pool Heights During a Decreasing Heat Flux 
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Figure 62. Performance Variation of Tubes 1 to 3 at a 10 cm 
Pool Height During a Decreasing Heat Flux with Auxiliary 
Heaters off (DTB06T) 
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Figure 63. Performance Variation of Tubes 1 to 4 at a 10 cm 
Pool Height During a Decreasing Heat Flux with Auxiliary 
Heaters off (DTB07T) 
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Figure 64. Performance of Tube 1 at a 20 cm Pool Height with the 
Auxiliary Heater Setting at 3kW 
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Figure 65. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 20 cm Pool Height 
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Figure 66. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 20 cm Pool Height 
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Figure 67. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 10 cm Pool Height 
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Figure 68. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 0 cm Pool Height 
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Figure 69. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 0 cm Pool Height 
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Figure 70. Nucleation Site Deactivation with a 1 kW Auxiliary 
Heater Setting and a 0 cm Pool Height 
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Figure 71. Nucleation Site Deactivation with a 3 kW Auxiliary 
Heater Setting and a 10 cm Pool Height 
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Figure 72. Nucleation Site Deactivation with a 3 kW Auxiliary 
Heater Setting and a 10 cm Pool Height 
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Figure 73. Nucleation Site Deactivation with a 3 kW Auxiliary 
Heater setting and a 10 cm Pool Height 
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Figure 74. Nucleation Site Deactivation with a 3 kW Auxiliary 
Heater Setting and a 10 cm Pool Height 
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Figure 75. Nucleation Site Deactivation with a 3 kW Auxiliary 
Heater Setting and a 10 cm Pool Height 
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Figure 76. Comparison of Tube 1 During a Decreasing Heat Flux for 
Smooth vs Turbo-B Bundles in R-113 
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Figure 77. Comparison of Tube 1 During a Decreasing Heat Flux for 
Turbo-B Tubes in a Bundle with R-113 vs a Bundle for R-114 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. When secured overnight, the bundle exhibited significant temperature overshoot 
for increasing heat flux (approx. 10°C). The hysteresis loop formed by increasing 
followed by decreasing heat flux was very repeatable. 

2. In the natural convection region, activating lower tubes within the bundle did not 
affect the upper tubes. However, the heat-transfer coefficient decreased as one moved 
down the bundle. 

3. Activation of the lower tubes caused the upper tubes to nucleate earlier (i.e. at a 
lower heat flux). As each tube successively nucleated, the heat-transfer performance of the 
remaining lower tubes improved. 

4. In the nucleate boiling region, the performance of the top tube was the best. This 
performance dropped off steadily as one moved down the bundle. 

5. Variation in liquid height above the bundle (from 0 to 20 cm) significantly delayed 
the point of incipience. 

6. For zero pool subcooling, nucleation site deactivation was greatly dependent on 
both pool height (i.e. pressure head) and tube turn off time. At zero pool height, 
deactivation was negligible for all time whereas for 10 and 20 cm, there was complete 
deactivation in 60 and 15 seconds respectively (i.e. decreasing with increasing pool 
height). 

7. At a typical heat flux of 20 kW/m^, the Turbo-B tube bundle heat-transfer 
coefficient was approximately four times that of a smooth tube bundle. 
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8. The uncertainty of the data was greatest at low heat flux and low wall superheat 
due to the measurement being of the order of accuracy of the thermocouples. 

B . RECOMMENDATIONS FOR FUTURE WORK 

1. Due to the way a real flooded evaporator operates, it would be preferable to heat 
the tubes within the bundle with hot water. 

2. The pool height tests should be repeated, but by varying the vapor pressure above 
the liquid, the local pressure at the tube surface should be maintained at a constant value. 

3. A pressure test should be conducted and a pressure relief valve fitted to enable the 
testing of alternative refrigerants such as R-124 which has a higher vapor pressure. 

4. Attention needs to be given to the question of refrigerant disposal. There are 
reclamation projects undertaken by most manufacturers; however, a method still needs to 
be found to remove the refrigerant from the apparatus into a container for reclamation. 

5. Neutrally buoyant particles should be placed in the pool to facilitate study of the 
circulation patterns within the bundle in more detail. 

6. A high speed camera should be used to study the nucleation process and 
circulation patterns in more detail (the mechanical engineering department is hoping to 
purchase such a camera in the near future). 

7. A sight glass should be fitted to the refrigerant storage tank. 
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Temp vs. Latent Heat 



APPENDIX B: 



PROPERTIES OF R-113 
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APPENDIX C: 



SAMPLE CALCULATIONS 



Data set number 1 of run DTB03 (Decreasing heat flux, Turbo-B tube. Surface 
preparation D) was used for the sample calculations in order to validate the program used 
for data acquisition DRP4A. The working fluid was R-1 13. 

1. Test tube dimensions 

= 11.60 mm 

Dq = 14.15 mm 
Dj = 12.70 mm 
L = 203.2 mm 
Lu = 25.4 mm 

2. Measured Parameters 
T1 =54.41 °C 

T2 = 55.25 °C 
T3 = 55.80 °C 
T4 = 55.15 °C 
T5 = 55.04°C 
T6 = 55.38°C 
Tldl=47.46®C 
Tld2 = 47.48®C 
Aas = 3.838V 
Vas = 3.958V 
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3. Calculations 



The heater power is first calculated from 

q=Vas(V) X Aas(V) x 60(VA^) x 1(AA0 

Note; The multiplication factors of volt and amp sensors are 60 and 1, respectively. 



Therefore: 



q= 3.958V(3.838V)(60VA^)(1A/V) 
= 911.45W 



The tube inside wall temperature is obtained from the average of all six thermocouple 
readings. 

_ 6 
Twi = ^ Tn 

6n=i 

= 1/6 (54.41 + 55.25 + 55.80 + 55.15 + 55.04 + 55.38) 

= 55.17°C 



The tube outside temperature is calculated by knowing the inside wall temperature and 
using Fouriers Conduction law . Uniform radial conduction is assumed. 



Two = Twi - 



q nn(^)l 

2Jt(kcu)(L) 
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where the second term on the right hand side is the Fourier conduction term. If we 



define this term as 
d> = 



q [ln(^)] 

2ti(kcu)^) 



and 



0b = Two - Tsatc 



where kcu is the thermal conductivity and is calculated as follow 
kcu = 434.0 - [0.1(Twi-K)] 

= 434.0 -[0.1(328.32)] 

= 401.17 W/mK 



Now, 



Two = Twi - 



911.45tln(|f^)l 

271(401.17X0.2032) 



= (55.17 - 0.3567)®C 



= 54.82®C 



The liquid saturation temperature at the top of the tube bundle is 

Tsat = 3ML±nd2. 

2 

^ 47.46+47.48 
2 

= 47.47°C 
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In order to calculate the local saturation temperature for each tube, correction factors 
are needed to account for hydrostatic pressure differences between the tube locations and 
the liquid free surface. This difference is calculated by: 

AP = p (g)(ht) 

where the density is calculated by 

p = 1620.7479 - Tsat(2.2 186346 + Tsat(2.3578291E-3)) 

Therefore, for Tsat = 47.47°C 
p= 1510.1 kg/m3 

So, for the top tube in the bundle which is 0.0124 m below the thermocouples 
measuring the pool temperature, 

AP = 1510.1 (9.81)(0.0124) 

AP = 183.69 Pa 

For a 183.69 Pa pressure difference,the corrected saturation temperature is obtained 
by adding 0.054°C (obtained from standard tables for R-113 or by equations shown in 
Appendix F.) to Tsat The corrected Tsat is: 

Tsatc = (Tsat + 0.054)°C 
= (47.47 + 0.054)°C 
Tsatc = 47.52°C 
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Therefore, the wall superheat can be obtained by the following 



6b = Two - Tsatc 
= (54.82 - 47.52)°C 
= 7.30°C 

Now that the wall superheat is known, we need to calculate the heat flux and the heat 
transfer coefficient. To do this, we know that the tube is 12 inches long and is heated in a 
eight inch center portion of the tube. The unheated lengths of the tube are a one inch and a 
three inch section on opposite ends of the tube. These unheated lengths have a fin effect 
during the heat transfer process to the evaporating refrigerant. In order to account for this, 
the following procedure was adopted for both the one and three inch sections. Calculations 
are shown below for the one inch section . 

Heat transfer from the unheated end is calculated as heat from the base of a fin: 

qr = (hb xp xkcu X Ac)°'^ x 0b x tanh(n xLc) 

where 

p = 7C X Do 
= 7 CX 14.15E-3 m 
= 44.45E-3m 

now 

Ac = J(d2-d?) 

Ac = ? (0.014152 - 0.012702) 

4 

Ac = 3.0578E-5 m^ 
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The corrected length of unenhanced surface at the end was calculated as follows 



Lc — Lu + 2 
U = 0.0254 

Lc = 0.0258 m 



hb is the natural convection heat transfer coefficient of the fin like ends and was 
calculated by using the Churchill-Chu [Ref. 24] correlation for natural convection from a 
smooth horizontal cylinder, as modified by Pulido [Ref. 27]. 






0.6 +0.387 X 



gxpxE^xObX tanh(n xLc) 



V X a X Lo X n 



l+r O.559 



Pr J 



16 



27 



where 



n = 



hbXp 



kcu xAc 



0.5 



Therefore an iterative technique was necessary to calculate hb. The iterative technique 

W 

used was to assume hb was 190 ^^2^ and to continue the iteration until successive values 
are within 0.001 of each other. The fluid physical properties are calculated at the vapor 
mean film temperature, given by the following equation. 
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rr Tsatc + Two 
Tfilm = 2 



Tfilm = 5I.17°C =324.32°K 



For R-1 13, the physical properties are given in the program by 



Dynamic viscosity, Tfiin, in ®K 

r 503 1 

|X = 1.34E-5 X lOLTfitoK - 2.15J 

|X = 4.88E-4 kg/m s 

Specific heat, Tfiin, ^ 

Cp = 929.0 + (1.03 X Tfilm) 

Cp = 981.71 J/kg-K 



Density, Tfiin, in ®C 

p = 1620.7479 - TfiU2.21 86346 + Tf,in,(0.0023578291)] 
p= 1501.1 kg/m3 

Thermal conductivity of R-1 13, Tfiin, ^ 
k = 0.071 - (0.000261 x Tfiim) 



k = 5.765E-2W/mK 



Prandtl number 



Pr = Cpxtt 



Pr = 981.71 



X 4.88E-4 
5.7645E-2 



Pr = 8.31 
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Thermal expansion 

p = -(l/p)(Ap/AT) 

P51.07 = 1501.29 kg/m3 Psi.27 = 1500.8 kg/m^ 

13 = -(1/1501.1X0.49/0.2) 

P= 1.63E-3 (1/K) 

Kinematic viscosity 

v = ^ 

P 

„ _ 4.88E-4 
1501.1 

V = 3.25E-7 m2/s 

Thermal diffusivity 

a = — ^ — 
pxCp 

g- 5.7645E-2 
1501.1 X981.71 

a=3.912E-8m2/s 

Knowing the above properties, the heat transfer coefficient, hb, can be obtained by 
iteration. 

hb = 185.42 W/m2.K 



knowing this we calculate n 



n = 



hbxp 
k{^ XAl<; 



0.5 
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n = 



185.42 X44.45E-3 
401.17 X3.0578E-5 

n = 25.92 

then we can obtain qf 

Qf = (hb X p xkcu X Ac)^‘^ X 6b X tanh(n xL^) 

qf = (185.42 X 44.45E-3 x 401.17 x 3.0578E-5)°-^ x 7.3 x tanh (25.92 x 0.0258) 
qf=1.36W 

The corresponding results for the three inch section are 
hb = 158.27 W/m2 K 
qf = 2.04W 

Therefore the heat transfer through the heated length of the tube is 
qs = q - qf(l inch section) - qf(3 inch section) 
qs = (911.45- 1.36-2.04) W 
qs = 908W 

and the heat flux and the heat transfer coefficient are as follows 



q" = qs/As 

= qs/(7ixDoXL) 

= 908/(71x0.01415x0.2032) 
= 908/9.033E-3 
= 1.005E+5 W/m2 
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and finally the heat transfer coefficient 



As X (Two-Tsat) 

= 908 / (9.033E-3 x 7.3) 
h = 1.376E+4 W/m2-K 
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APPENDIX D; UNCERTAINTY ANALYSIS 



The same data run (DTB03) was chosen for the uncertainty analysis. Therefore, the 
measured and calculated parameters found in the sample calculation were used in this 
section. The uncertainty analysis performed was for a high heat flux but the same 
procedure could be performed at any heat flux to determine the uncertainty bands. All 
uncertainties are presented as a percentage of the calculated parameter. The uncertainty 
associated with the experimental parameters is calculated from the equation suggested by 
Kline and McClintock [Ref 28]. For example, if 
R R • • • » 

then 

5R = 

where: 

5R = uncertainty of the desired dependant variables 
Xn = measured variables 
5xn = uncertainty in measured variables 



,3R- / ,3R_ / 3R_ 

(t^xi ) + (:r^x2 ) + • • • + (:~^xn ) 
OXi 0X2 OXn 



0.5 



The boiling heat transfer coefficient is given by 



h = 



Qs 



As X (Two-Tsat) 



where 



_ _ q [Intel 

Two = Twi- „ tte 
27i(kcu)(L) 
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In the above equation, the second term on the right hand side is usually called the 



Fourier heat transfer conduction term. If we define this as O, then 

q [ln(fel 

Dtp - 

2;t(kcu)(L) 



and 



6b = Two - Tsatc 



With this notation, the uncertainty in the heat transfer coefficient is obtained using the 
following equation. 



5h = 



h 



]5q' 


2 

+ 


’6As‘ 


2 

+ 


STwi 


2 

+ 




2 

+ 


§Tsat 


2 


_[q. 




[AsJ 




9b J 




.9b. 




9b J 





where 



q = V X I 

q = V(V) X I(V) X 60(VAO X 1(AA0 



and the uncertainty is 



§ 9 .= 

q 



5Y 
li VJ 




1 

2 



The accuracy in the voltage and current sensors are as follows 



5Vas = ±0.05 V 6Aas = ±0.025 A 
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From the sample calculation section 

Vas = 3.958V Aas = 3.838V 

Therefore, 

5q/q = ((SVasA^as; + (5Aas/Aasr) 

5q/q = ((0.05/3.958)2 + (0.025/3.838)2) 0.5 
5q/q = 1.42 percent 



Calculation of the surface area and the uncertainty of it are done as follows 
As = JUxDoXL 



5As 

As 









Knowing the dimensions of the tube from the manufacturer and estimated 
inaccuracies from work shop tools and human error, the uncertainty was calculated. 



Dimensions 

Do= 14.15 mm L = 203.2 mm 



Inaccuracies in measurements 

5Do = 0.1 mm 5L = 0.2 mm 
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Uncertainty analysis performed 

o o 0-5 

5As/As = ((5Do/Dor + (5L/Lr> 

5As/As = ((0.1/14.15)2 + (0.2/203.2)2)°-^ 
5 As/As = 0.7135 percent 



The uncertainty calculation for the Fourier conduction term is given below 

1 

2 



M = 




2 


5 kcu] ^ 


2 


<D 


..q. 




L kcu J L L J 


- 



kcu was calculated using 

kcu = 434.0 - [0.1(Twi-K)l 



= 434.0- [0.1(328.32)1 
= 401.17 W/mK 



and its uncertainty 

5kcu = [(0. 1 X 5Twi-K)^] 2 



6Twi and 5Tsat are obtained using uncertainties in the thermocouple readings. 
Average wall inside temperature Twi was obtained taking the average of six thermocouple 
readings inside the tube. The uncertainty associated with this variable is 



5Twi = 





where 5Tc was obtained as ± 0.5®C from "Omega 1987 Complete Temperature 
Measurements Handbook pp. T-37. Using engineering judgement this seems high will use 



o.rc. 
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5Twi = ^6 ^ 



5Twi = 0.04 °C 



Saturation temperature was obtained by taking the average of two thermocouple 
readings and the uncertainty in this temperature was calculated from the following equation. 



5Tsat = 


[2x1 


If] 


11 


5Tsat - 


[2x| 




1 



5Tsat = 0.07®C 



Knowing the uncertainty in the temperatures, we can now calculate the uncertainties 
in the following. 

5kcu=[(0.1 x5Twi-K)^]2 
6kcu=[(0.1 X 273.19)^] 2 
6kcu = 27.32 W/mK 



Now we can calculate the uncertainty in the Fourier conduction term 



SO = 


6q 


2 

+ 


Skcu 


2 


sl' 


2 


<D 


..q. 




. kcu - 




LlJ 


. 



= [[0.0142] 
d) ^ 



21 



401 



22_1%[JL2. 

M\ l203.: 
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= 6.96 percent 

from the sample calculations we know that <I> = 0.3567 °C or 
5<D = 0.3567 X 6.96 percent = 0.025 



therefore; 



]6q' 


2 

+ 


‘5As' 


2 

+ 


5Twi 


2 

+ 


'50 


2 

+ 


5Tsat 


2 


[q. 




. As . 




0b J 




.6b. 




0b J 





^ = [[0.0142] 2 + [0.07 1 35] U [^] ^ + [0^] ^ + [!^] 2 

^ = 7.4 percent 
h 



Table 3 shows the results of the uncertainty analysis performed. The high and low 
heat flux correspond to the approximate values of 1.005E+5 W/m^ and 597 W/m^ 
respectively. 
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TABLE 3. 



UNCERTAINTY ANALYSIS RESULTS 



Variable 


High Heat Flux 


Low Heat Flux 


AT 


7.3 


0.38 


Twi 


55.17 


48.04 


Tsat 


47.47 


47.6 


5V 

V 


1.26% 


16.4% 


§l 

I 


0.65% 


8.3% 


5q 

q 


1.42% 


18.% 


5Do 

Do 


0.707% 


0.707% 


5L 

L 


0.098% 


0.098% 


5As 

As 


0.714% 


0.714% 


5kcu 

kcu 


6.8% 


6.8% 


<D 


6.96% 


21.9% 


5h 

h 


7.4% 


22.4% 
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APPENDIX E: OPERATING PROCEDURES 



SYSTEM MGHT-OFF 

1 . Power to the 8 ton refrigeration unit is provided by the breakers located in 
the main distribution panel located in the laboratory. These breakers were never secured. 
However, if power to this panel was lost then these breakers must be reset. 

2. Turn the switch on the refrigeration unit control panel, located in front of the 
refrigeration unit to the "auto" position after passing through "on" position. This switch is 
always left on, unless unit was taken down for long repairs. 

2. Push the start button in the control box for the recirculation pump. This 
control box is located on the bulkhead above the recirculation pump in the outside area 
adjacent to the refrigeration unit 

3. Set the desired temperature on the roughly graduated Fahrenheit scale on the 
control panel thermostat. It requires approximately one hour to chill the sump 5°C. The 
thermometer of the ethylene glycol sump must be monitored to ensure the desired sump 
temperature is attained and maintained. Slight adjustments in the refrigeration unit 
thermostat can be expected due to the coarseness of its scale. 

4. Energize the desired pumps by switching on the breakers in the main 
distribution power panel. Once the power is energized in the main power panel turn on the 
pump motors by pushing down on the arm of the appropriate breaker box for the pumps 
located on the bulkhead next to the ethylene glycol sump. The pumps are marked 
"auxiliary condenser"(Pump #2) and "instrumented tube condenser"(Pump #1), 
respectively in the breaker box. Flow in the auxiliary condensate system can be controlled 
with the individual globe valves located at the coil penetrations on the apparatus. The 
auxiliary condenser will produce the fastest adjustments to system pressure if pressure is 
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rising too quickly. Row through the instrumented tubes can be controlled by the ball 
valves located at the bottom of the respective flow meter. 

5. Energize the heater variacs desired by switching on the breakers (Auxiliary, 
Simulation, and/or Bundle) in the main distribution power panel (located near entrance of 
lab) and individual breakers for each of these in the power distribution box (near 
apparatus). 

6. After ensuring that the breakers for the heater tubes desired are in the "on" 
position, adjust variac/s to the desired amount of heat flux. Monitor apparatus pressure 
through the pressure gage at the top of the apparatus, ensuring system pressure does not 
exceed 30 psig. 

SYSTEM SHUTDOWN 

1 . Turn all variacs to the zero position and switch off all breakers in the power 

panels. 

2. If apparatus will not be operated for an extended period, turn the switch on 
the refrigeration control panel to the "off position after passing through "on". 

3. Allow the recirculation pump to operate for at least five minutes after 
switching off the refrigeration unit to dissipate any back pressure in the system. 

4. Turn the breakers for the pumps to the off position at the switch boxes near 
the ethylene glycol sump, and then secure power at the main distribution power panel. 

EMERGENCY SHUTDOWN 

1. Secure all power at the main distribution power panel 

2. Pull fire alarm 

3. Evacuate building 
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APPENDIX F: 



PROGRAM DRP4A 



10001 
1004 I 
1008« 
10121 
101G 
1020 
1024 
1028' 
1032 
103G 
1040 
1044 
1048 
1052 
105G 
10G0 
10G4 
10G8I 
10721 
107G 
1080 
1084 
1088 
1092 
109G 
1100 
1 104 
1108 
1112 
1116' 
1120 
1124» 
1128 
1132 
Ldtc( 
1136 
11401 
11441 
1148 
1152 
115G 
11G0! 
1 1G4 I 
11G8 
1172 
1 17G 
11801 
1184 
1188 
1192 
119G 
1200 
1204 
1208 
1212 
121G ' 



FILE NAME^ DRP4A 

DATE- November 22, 1988 

REUISED: May 1991 (S. MEMORY) 



BEEP 

PRINTER IS 1 
I dp=0 

PRINT USING "4X, *• 
PRINT USING "GX, " 
PRINT USING “GX, " 
PRINT USING "GX , " 
PRINT USING "GX, " 
PRINT USING "GX," 
PRINT USING "GX,"' 
PRINT USING "GX," 
PRINT USING "GX,"' 



Select option default is 0-""" 

0 Taking data or re-process i ng previous data" 

1 Plotting data on Log-Log 

2 Plotting data on Linear 

Purge " 



Fixup 

Move 

Comb 

Read Plot""" 



IDP IS A PROGRAM VARIABLE TO SELECT A SUBROUTINE 
INPUT Idp 

IF ldp=0 THEN CALL Main 
IF Idp=l THEN CALL Plot 
IF Idp*2 THEN CALL Plin 
IF Idp»3 THEN CALL Purg 
IF Idp*4 THEN CALL Fixup 
IF Idp=5 THEN CALL Move 
IF Idp*G THEN CALL Comb 
IF Idp=7 THEN CALL Readplot 
END 

SUB Mam 

ICAL=THERMOCOUPLE CALIBRATION 
COM /Cc/ C(7> 

DIM Emf(35),T(35 >,Dla(G ),D2a(G ) , Dia( 6 ) , Doa( G > , La( G > , Lua( G > , Kcua( G ) ,Et ( 19 > , 
!0 >,Volt(2 ),Amp( 11 >,Twa(5 >,Tw(5 >,Theta(5),Thetab(5 ),0(5 >.01(5 >,0dp(5 ) 

DIM Htube(5> 

THERMOCOUPLE ARRAY (C( )) INITIALIZATION 

DATA 0. 1008G091 ,25727.94569, -767345.8295,79025595. 81 
DATA -9247486589,6 .97G88E+ 11,-2. GG192E+ 13,3. 94078E+ 14 
READ C( * > 

PRINT HEADER AND INITIALIZE TIME CLOCK 

PRINTER IS 701 

BEEP 

INPUT "ENTER MONTH, DATE AND TIME ^MM^DD HH^MM 5S)",Date$ 

OUTPUT DIRECTED TO DATA A0UI5ITI0N SYSTEM (HP 3497A ) 

OUTPUT 709: "TD" ;Date$ 

OUTPUT 709; "TD" 

ENTER 709;Date$ 

PRINT 

PRINT " Month, date and tine "iDateS 

PRINT 

PRINT USING "10X,""NOTE Program name • DRP4 

BEEP 
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1Z20I DN IS THE UARInBLE FOR DISC NUMBER FOR RECORD KEEPING ONLY 

1224 input "ENTER DISK NUMBER", Dn 

1228 PRINT USING "lBX,""Disk number * "“,ZZ";Dn 

1232 BEEP 

1236 lm=0 

1240 INPUT "ENTER INPUT MODE ( 0=3497A , 1-FILE ) 0-DEFAULT" , Im 
1244! 

1248! INPUT MODE ZERO IS FROM THE DATA AQUISITION SYSTEM 
1252 IF lm*0 THEN 
1256 BEEP 

1260 INPUT "GIUE A NAME FOR THE RAW DATA FILE",D2f i le$ 

1264 PRINT USING ** 16X , “ "Fi le name-- " " , 14A" ; D2f i le$ 

1268' 

1272! CREATE BDAT FILE ON THE MASS STORAGE MEDIA 

1276 CREATE 8DAT D2file$,B0 

1280! CREATE AN INPUT/OUTPUT LINK TO OPEN FILES 

1284 ASSIGN @File2 TO D2file$ 

1288J 

1292! CREATE DUMMY FILE UNTIL Nrun KNOUN 

1296 Dlfile$-"DUMMY- 

1300 CREATE 8DAT Dlfile$,60 

1304 ASSIGN @Filel TO DlfileS 

1308 OUTPUT @FileliDate$ 

1312! 

13161 CREATE A PLOT FILE 

1320 8EEP 

1324 INPUT "GIUE A NAME FOR THE PLOT FILE",Pfile$ 

1328 CREATE BDAT Pfile$,30 

1332 ASSIGN @Plot TO PfileS 

1336 8EEP 

13401 

13441 IDTC - NUM8ER (TOTAL) OF DEFECTIUE THERMOCOUPLES 

1348 INPUT "ENTER NUMBER OF DEFECTIVE TCS ( 0-DEFAULT ) Mdtc 

1352! LDTC - LOCATION OF DEFECTIVE THERMOCOUPLE 

1356* 

1360 IF Idtc-0 THEN 

1364 PRINT USING "l6X,“"No defective TCs exist""’ 

1368 ELSE 

1372 PRINT USING " 16X "Defect ive Thermocouples Indicated by -99.99""" 

1376 END IF 

1380! 

1384 BEEP 

1388! DEFECTIVE THERMOCOUPLES MAY 8E IN CHANNELS 40-69 

1392! THEMOCOUPLES ARE ENTERED AS DEFECTIVE 8Y COMPUTER CHANNEL NO. 

1396! JDTC-COUNTER IN LOOP FOR DEFECTIVE THERMOCOUPLES 

1400* 

1404 IF ldtc>0 THEN 

1408 FOR Jdtc=0 TO Idtc-1 

1412 INPUT "ENTER DEFECTIVE TC LOCATION (BY COMPUTER CHANNEL NUMBER) 

" ,Ldtc( Jdtc ) 

1416 BEEP 

1420 NEXT Jdtc 

1424 END IF 

1428 PRINTER IS 701 

1432 OUTPUT @Filel ;Ldtc( O 

1436 * ^ 

1440 ! Im=l option (THIS OPTION ALLOWS DATA ENTRY WITH DATA FILE) 

1444 ELSE 
1448 BEEP 

1452 INPUT "GIVE THE NAME OF THE EXISTING DATA F ILE " , D2 f i 1 e$ 
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1456 

1460 

1464 

1468 

1472 

1476 

1480 

1484 

1488 

1492 

1496 

1500 

1504 

1508 

1512 

1516 

1520 

1524 

1528 

1532 

1536! 

1540 

1544 

1548 

1552 

1556 

1560 

1564 

1568 

1572 

1576! 

1580 

1584 

1588 

1592 

1596' 

1600 

1604 1 

1608 

1612 

1616 

1620 

1624 

1628! 

1632 

1636! 

1640 

1644 

1648 

1652 

1656 

16601 

1664 I 

1668 

1672 

1676 

1680 

1684 

1688' 

1692' 



PRINT USING *’ 16X , ""Fi le name " 'M 4A " ; D2 f i 1 
ASSIGN ®File2 TO D2file$ 

ENTER @File2:Nrun 

ENTER @File2;Dold$,Ldtc( * ), Itt ,Bop,Nht,Natp,Nrt , Corn 
BEEP 

INPUT "GIUE A NAME FOR PLOT FILE\Pfile$ 

CREATE 8DAT Pfile$,30 
ASSIGN ©Plot TO PfileS 

PRINT USING "lGX.""This data set taken on ■ " ’M 4A " . Dol d$ 

8EEP 

PRINTER IS 1 

PRINT USING MX /'"SELECT TU8E TYPE""" 

PRINT USING "6X/"0 SMOOTH 

PRINT USING "6X,""1 FI NNED( 1 9/ IN > 

PRINT USING "6X/'“2 HIGH FLUX """ 

PRINT USING "6X,""3 TUR80-B """ 

INPUT Itt 
END IF 

IF In=l THEN GOTO 1768 
PRINTER IS 1 

IF ln=0 THEN 



PRINT 


USING “4X/‘ 


"Select tube type""" 


PRINT 


USING "EXC- 


" 0 


Smooth 


PRINT 


USING "EXC- 


" 1 


FINNED 19/IN (DEFAULT)" 


PRINT 


USING "EXC- 


" 2 


HIGH FLUX""" 


PRINT 


USING "6Xc“ 


- 3 


TURBO-B--" 


PRINT 


USING "EXC- 


" 4 


GROUTH""" 


PRINT 


USING "6Xc" 


- 5 


GROWTH""" 


PRINT 


USING "6Xc" 


" 6 


GROWTH 


ITT=TUBE TYPE 






INPUT 


Itt 






OUTPUT @Filel;Itt 







END IF 

PRINTER IS 701 
Itt-2 

PRINT USING "lBX,""Tube Type; "",DD";Itt 

BEEP 

8op=0 

INPUT "ENTER 8ULK OIL % (DEFAULT-0%> ",Bop 
OUTPUT @Filel;8op 

PRINT USING "16X,""Bulk oi IX*" " , DD" ; Bop 
8EEP 

NHT=NUM8ER OF ^ HEATED TUBES 
Nht=5 

INPUT "Enter number of heated instrumented tubes( de f aul t =5 ) " , Nht 
OUTPUT ©Filel ;Nht 

PRINT USING " 16X “Number of heated instrumented tubes* "" ,DD" ; Nht 
8EEP 



Natp*Number of active dummy pairs 
Na t p*0 

INPUT "Enter number of active dummy pairs ( De f au 1 t *0 > " , Na t p 
OL'TPUT ©Filel;Natp 

PRINT USING " 1 BX /'" Number of active dummy pairs*" ", DD" . Natp 
8EEP 

NPT-NUM8ER OF ADDED HEATED TU8ES TO ENHANCE BUNDLE EFFECT 
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1696 Nrt=0 

1700 INPUT "Enter number of added heated tubes from simulation heat ers ( Def aul t » 
0>" ,Nrt 

1704 OUTPUT @Filel;Nrt 

1708 PRINT USING " 1 6X Number of added heated tubes( from simulation heaters)*" 
■■ ,DD" ;Nrt 
1712 BEEP 
1716! 

17201 CORR IS CORRECTION FOR INSTRUMENTED TUBE HEIGHT 
1724 Corr*0 

1728 INPUT "WANT TO CORRECT TSAT FOR TUBE HEIGHT ( 0*YES( DEFAULT ), 1*N0 >", Corr 
1732 IF Corr*0 THEN PRINT USING "16X,""TSAT is corrected instrumented heal 

ed tube height 

1736 IF Corr*l THEN PRINT USING "16X,"“TSAT is NOT corrected for instrumen 

ted heated tube height""" 

1740 OUTPUT @Filel;Corr 
1744 8EEP 

1748! ILQU-INPUT MODE^ LIQUID, UAP0R,0R LIQUID UAPOR AUERAGE 
1752 Ilqv-0 

1756 INPUT "SELECT ( 0*LIQ( def aul t ) , 1 -UAP , 2 = ( LIQ+UAP >/2 ) " , I Iqv 
1760! 

1764! DlA*Diameter at thermocouple positions (meters) 

1768 DATA .0122,0.0098,0.0106,0.0116,0,0,0 
1772 READ Dla(*) 

1776 Dl=Dla(Itt) 

1780! 

1784! D2*Diameter to base of fins (outside dia for smooth )( met ers ) 

1788 DATA .0158,0.0125,0.0158,0.01415,0,0,0 
1792 READ D2a( * > 

1796 D2»D2a(Itt) 

1800! 

1804! Di*Inside diameter of unenhanced ends (meters) 

1808 DATA .0132,0.0109,0.0116,0.0127,0,0,0 
1812 READDia(«) 

1816 Di-Dia(Itt) 

1820! 

1824! Do*0utside diameter of unenhanced ends (meters) 

1828 DATA .015875,0.0125,0.015875,0.01415,0,0,0 
1832 READ DoaC) 

1836 Do*Doa(Itt) 

1840! 

1844! L“Length of enhanced surface (meters) 

1848 DATA .2032, .2032, .2032, .2032, .2032, .2032, .2032 
1852 READ La( * ) 

1856 L=La(Itt) 

18601 

1864J Lu*CORRECTED Length of unenhanced surface at the ends (METERS) 

1868ILU=LFIN + THICKNESS/2 

1872 DATA . 0261 0254 ,. 0264 , 0 . 0258 , 0 , 0 , 0 

1876 READ Lua( * ) 

1880 Lu=Lua( I tt ) 

18841 

1888! LU=corrected length of 3 inch finned like end 
1892 DIM Lva(6) 

1896 DATA . 0769 ,. 0762 ,. 0773 , 0 . 0766 , 0 , 0 , 0 
1900 Rf AD Lva( * ) 

1 904 Lv*Lva( 1 1 1 ) 

1908I Kcua=Thermal Conductivity of tube 
1912! DATA 401,0,0,0,0,0,0 
1916' READ Kcua( ♦ ) 
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1920» 


Kcu=Kcua( 1 1 1 ) 




1924 


A=PI*(Do"2-Di '2 )/4 


1928 


P=PI*Do 




1952 


J = 1 




1956 


5x=0 




1940 


Sy=0 




1944 


Sxs=0 




1948 


Sxy=0 




1952 


Repeat- ' 




19561 






1960 


IF In-0 THEN 




1964* 


Ot 1 d*des ired 


temperature of liquid 


1968 


Dtld=47.5 


IR-113 


19691 


Dtld=2.2 


‘R-114 


1972 


I do-2 




1976 


ON KEY 0,15 


RECOUER 1952 


1980 


PRINTER IS 1 




1984 


PRINT USING 


"4X, ""SELECT OPTION """ 


1988 


PRINT USING 


"6X,""0-TAKE DATA""" 


1992 


PRINT USING 


"6X, ""1-SET HEAT FLUX 


1996 


PRINT USING 


"6X, ""2-SET Tsat (DEFAULT SET 


2000 


PRINT USING 


"4X,""N0TE: KEY 0 = ESCAPE""" 


2004 1 


I do=des ired 


opt ion 


2008 


BEEP 




2012 


INPUT Ido 




2016» 






2020 


BEEP 




2024 1 


Set default 


value for input 


2028 


IF Ido>2 THEN Ido=2 


2032! 


Take data option 


2036 


IF Ido-0 THEN 2440 


2040! 






2044 1 


LOOP TO SET HEAT 


FLUX (FOR TOP INSTRUMENTED 


2048 


IF Ido-1 THEN 


2052 


Dqdp-100000 





2056 
Otot • 
2060 
(W)**' 
2064 
20681 
2072 ‘ 



PRINT USING “4X, ""Qdp 



PRINT USING "4X/’*‘(U/fn"2) 



QDPs in 



( W/n"2 > 



Nrt 



Err*l 

Reset, read channel 25-30, autonat ic scaling 

Channel 25*aux anp5,2G=sin anps,27*inst volts, 28*sin 



vol ts , 30=inst anps 



2076 

2080 

2084 

2088 

2032 

2096 

2100 

2104 

2108 

2112 

2116 

21201 



OUTPUT 709; "AR AF25 AL50 UR5" 
FOR 1=10 TO 11 

OUTPUT 709; "AS SA" 

ENTER 709; Amp( I ) 

NEXT I 

FOR 1=0 TO 2 

OUTPUT 709. "AS 5A" 

ENTER 709;Uol t( I ) 

NEXT I 

OUTPUT 709: "AS SA" 

ENTER 709.Anp(0) 

Calculate actual heat flux 



2124 , O(0>=60*Uolt(0)*Anp(0) 

2128 Qdp(0)=Q(0)/(PI*O2*L) 

2132 Osin=60*20*Uolt( 1 )*Anp( 11 ) 

2156 Qdp5ifn=0sifn/( PI *02 *. 2032*5 ) 

2140 Oau>.=60*20*Uol t ( 2 >*Anp( 10 ) 



Qdpaux 
(U/n''2 ) 



volts , 29*aux 
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2144 

2148 

2152 

2156 

2160 

2164 

2168 

2172 

2176 

2180 

2184 

2188 

2192 

( 0 ) fOdps 

2196 

2200 

2204 

2208 

2212 » 

2216! 

2220 

2224 

2228 

2232 

2233' 

2236 

2240 

2244 

2248 

2252 

2256 

2260 

2264 

Tvat 

2268 
2272! 
2276 
22801 
rnocoupl 
2284 
2288 
2292 
2296 
2300 
2304 
2308 
2312 
2316 
2320 
2324 1 
2328 
23321 
2333 



i Odp 



(DEFAULT-47.5 C - R“ll3)“,Dtld 
(DEFAULT-2.2 C - R-114)\Dtld 



Odpaux*Qaux/(PI*.0160*. 1776*4 ) 

0tot=0( 0 )*Nht +0s in+Oaux 
Nrt-0dpsin/0dp( 0 ) 

IF ABS( Aqdp-Dqdp )>Err THEN 
IF Aqdp>Dqdp THEN 
BEEP 4000,. 2 

ELSE 

BEEP 250,. 2 
END IF 

IF Nrt<.l THEN Nrt-0 
IF Odpaux<100 THEN Qdpaux-0 
IF Odp5in<100 THEN Odp5iPi-0 

PRINT USING “4X,2(MZ.3DE,2X),2X,(MDD.DD),2X,2(MZ.3DE,2X)“ 
xn , Nrt , Qdpaux , Otot 
UAIT 2 
GOTO 2076 
END IF 
END IF 



LOOP TO SET Tsat 
IF Ido*2 THEN 

IF Ikdt-1 THEN 2240 
BEEP 

INPUT "ENTER DESIRED Tsat 
INPUT -ENTER DESIRED Tsat 
Ikdt«l 
Oldl»0 
Old2-0 
Nn-1 

Nrs-Nn MOD 15 
Nn-Nn+ 1 
IF Nrs-l THEN 

PRINT USING -4X,“- DTsat Tldl Tld2 Tlbb 

Tvab Tlav 

END IF 

Read thermocouple voltages for vapor, liquid 
OUTPUT 709; "AR AF0 AL5 UR5" 

Sample each thermocouple 20 times and report temp for each the 
, vapor-0,1,2; liquid-3&4 
FOR 1-0 TO 5 
Sum-0 

OUTPUT 709; -AS SA- 
POR Ji-1 TO 20 

ENTER 709;Eliq 
Sum-Sum+El iq 
NEXT Ji 
Emf ( I )=Sum/20 
T( I )-FNTvsv(Emf( I )) 

NEXT I 

Compute average temperature of liquid 
Tlav=(T(3) + T(4 ))*.5 

Compute average temperature of vapor 
Tvavl-(T(0)+T( 1 ) )/2 



2334 

2336 

2340 

2344 

2348 

2552 

2356 



Tvav2=T(2 ) 

Tvav = (T(0 )+T( 1 ) + T( 2 ) )/3 
IF ABS(Tlav-Dtld)>.2 THEN 
IF Tlav>Dtld THEN 
BEEP 4000, .2 

ELSE 

BEEP 250,. 2 
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END IF 



2360 

2364 

2368 

2372 

2376 

2380 

2384 

2388 

2392 

2396 

2400 

2404 

2408 

2412 

2416 

V 

2420 

2424 

2428 

24321 

2436' 

2440 

2444 

2448 

2452 

2456 

2460 

2464 

2468 

2472 

2476 

2480 

2484 

2488 

2492 

2496 

2500 

2504 

2508 

2512 

2516 

2520 

2524 

252B 

2532 

2536 

2540 

2544 

2548 

pea t 

2552 

2556 

2560 

2564 

2568 

2572 

2576 

2580 

2584 

2588 



ELSE 

IF ABS(Tlav-Dtld)>.l THEN 
IF Atld>Dtld THEN 
BEEP 3000, .2 

ELSE 

BEEP 800, .2 
END IF 
END IF 
END IF 

Errl-Tlav-Oldl 
Oldl-Tlav 
Err2=Tvav-0 1 d2 
01 d2=»Tvav 

PRINT USING "4X,7(MDDD.DD,3X )";Dtid,T(3),T(4 ) , T( 5 ) , Tvavl . Tvav2 , Tla 

WAIT 2 
GOTO 2252 
END IF 

TAKE DATA IF ln=0 LOOP 
IF Ikoi=l THEN 2452 
BEEP 
Ikoi=l 

OUTPUT 709; "AR AF0 AL5 UR5" 

FOR 1=0 TO 5 

OUTPUT 709; "AS SA" 

Sun=0 

FOR Ji=l TO 20 
ENTER 709; E 
Sum*Sun+E 

IF I>2 THEN Et( Ji-1 )«E 
NEXT Ji 
Kdl-0 

IF I>2 THEN 

Eave=Sun/20 
Sun=0 . 

FOR Jk*0 TO 19 

IF ABS<Et( Jk )-EaveX5.0E-6 THEN 
Sun*Sun+Et ( Jk ) 

ELSE 

Kdl=Kdl+l 
END IF 
NEXT Jk 

IF I>2 THEN PRINT USING "4X,""Kdl = "'\DD";Kdl 
IF KdlM0 THEN 
BEEP 
BEEP 

PRINT USING ‘'4X,""Too nuch scattering in data - re 

data set 

GOTO 1980 
END IF 
END IF 

Enf ( I ) = Sun/( 20-Kdl ) 

NEXT I 

, OUTPUT 709; "AR AF40 AL69 UR5" 

FOR 1=6 TO 35 

OUTPUT 709; "AS SA" 

Sun=0 

FOR Ji=l TO 5 
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2592 




ENTER ■^09;E 


2596 




Sur=Surn + E 


2G00 




NEXT Ji 


2604 




Enf( I ) = Sun/5 


2608 




NEXT I 


26121 






2616! 




READ VOLTAGES ( 27=Ins t , 28*Sin, 29*Au: 


2620 




OUTPUT 709; "AR AF27 AL29 UR5“ 


2624 




FOR 1=0 TO 2 


2628 




OUTPUT 709; "AS SA" 


2632 




ENTER 709;Uolt( I ) 


2636 




NEXT I 


2640! 






2644! 




READ CURRENTS (30-34=In5t tubes; 35-: 


2648 




OUTPUT 709; "AR AF30 AL39 UR5" 


2652 




FOR 1=0 TO 9 


2656 




OUTPUT 709; “AS SA" 


2660 




ENTER 709; Anp( I ) 


2664 




NEXT I 


26681 


Read Current s( 25=Aux anp5,2G*Sin amps) 


2672 




OUTPUT 709; "AR AF25 AL2G UR5" 


2676 




FOR 1=10 TO 11 


2680 




OUTPUT 709; "AS SA“ 


2684 




ENTER 709;Anp( I ) 


2688 




NEXT I 


2692 




ELSE 


2696 




ENTER @File2 ;Enf< * ) .Uolt ( * ) , Anp( * ) 


2700 


END 


IF 


2704 1 






2708' 


CONVERT EMF’S TO TEMP , VOLT .CURRENT 


2712 


FOR 


1=0 TO 35 


2716 




T( I )=FNTvsv(Enf( I ) ) 


2720 




IF I>5 AND ldtc>0 THEN 


2724 




FOR li*0 TO Idtc-1 


2728 




IF Ldtcdi ) = I-5+39 THEN T( I 


2732 




NEXT Ii 


2736 




END IF 


2740 


NEXT I 


27441 


Nt c = 


*nr of thermocouples 


2748 


Nt c = 


*6 


2752 


FOR 


li*0 TO 4 


2756 




Q( Ii )=G0*Uolt(0)»Amp(Ii ) 


27601 




Tua*Averaoe temperature of the wall 


2764 




Tua( Ii )*0 


2768 




Ndtc=0 


2772 




FOR 1=1 TO Ntc 


2776 » 




Nn is counter in temp array, st( 


ernocouple in 


the tube ban^ ) 


2780 




Nn= I 1 *6+ 1+5 


2784 




IF ABS(T(Nn)) 99 THEN 


2788 




T(Nn >=-99.99 


2792 




Ndtc=Ndt c+1 


2796 




ELSE 


2800 




Twa( Ii ) = Twa( Ii ; + T(Nn) 


2804 




END IF 


2808 




NEXT I 



2812 Tuia( Ii ) = Twa< Ii )/(6-Ndtc ) 

281G NEXT Ii 

2820 Tlav=(T(3>+T(4))/2 

2321 Tvav=T(2) 
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:s:4 1 
2328 » 
2829' 
2832 
2836 
28401 
2844 I 
2848 
2852 
2856 
2860 
2864 
2868 
2872' 
2876 
2880 
2884 1 
2888! 
2892 
2896 
2900 
2904 
2908 
2912 
2916 
2920 
2924 
2928! 
29521 
2936 1 
2940 
2944 
2948 
2952 
2956 
2960 
2964 
2968 
2972 
2976 
2980 
2984 
29881 
2992 * 
2996 
3000 
3004 
3008 
3012 
3016 
3020 
3024 
3028 
5032 
3036 
3040 
3044 
3048 
3052 
3056 



Tvav = (T^0;**'T( 1 ;+rt2)>/5 



Tlav=T(5) 

Tcu = Tua( 0 ) 

Kcu=FNKcu( Tcu ) ' THERMAL CONDUCTIVITY OF COPPER 

I IF CURVE FIT NOT AVAIL USE ARRAY KCU( * ) 
FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
FOR 1=0 TO 4 

Tu( I )*Tua( I )-Q( I )*LOG(D2/Dl > / ( 2 *P I *Kcu *L ) 

IF llqv=0 THEN Texs=Tlav 

IF Ilqv=l THEN Texs=Tvav 

IF Ilqv=2 THEN Texs = ( T lav + T( 2 ) )* . 5 
IF Corr=l THEN Thetab( I )=Tw( I )-Texs 

IF Corr=0 THEN Thet ab(I ) = Tu( I )-( Tca s+ . 056+1 * . 1 29 > iR-114 

IF Corr=0 THEN Thet ab( I ) = Tw(I )-( Tex5+ . 054 + 1 ♦ . 144 ) IR-113 

NEXT I 



COMPUTE VARIOUS PROPERTIES 



Tfilm*=(Tu(0> + Tex5 )♦. 
Rho=FNRho( Tf i In ) 
Mu=FNMu ( Tf 1 1 m ) 
K=FNK(Tf iln) 

Cp = FNCp( Tf 1 In ) 
8eta=FN8eta( Tfiln ) 

N 1 =Mu/Rho 
Alpha*K/( Rho*Cp ) 
Pr*=Ni /Alpha 



5 ‘FILM TEMPERATURE 
'DENSITY 
•VISCOSITY 

•THERMAL CONDUCTIVITY 
•SPECIFIC HEAT 
•THERMAL EXPANSION 
IKINEMATIC VISCOSITY 
•THERMAL DIFFUSIVITY 
•PRANDTL 



COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
FOR UNENHANCED END( S > 

Lu=Lua( Itt > 

Hbar*190 

Fe = (Hbar*P/(Kcu*A ) >''.5*Lu 
Tanh = FNTanh( Fe > 

Theta(0)=Thetab(0)*Tanh/Fe 

Xx*=( g.81*Beta*Thetab( 0 )*Do''3*Tanh/( Fe*N i*Alpha ) )". 166667 
Yy=( l + ( .559/Pr )"( 9/16 ) >"(8/27 ) 

Hbarc=K/Do*( . 6+ . 387*Xx/ Yy )"2 
IF A8S( (Hbar-Hbarc )/Hbarc >>.001 THEN 
Hbar=( Hbar+Hbarc )*.5 
GOTO 2948 
END IF 



COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) 
QK0 > = ( Thetab(0 )*Tanh )*( ( Hbar*P*Kcu*A )'‘ .5 ) 
Oq=Ol(0 )+Oq 
Z = Z + 1 

IF Z=1 THEN 
Lu=Lv 
GOTO 2944 
END IF 
Z=0 

Olpct =Oq/Q( 0 ) 

Oq=0 

As=PI*D2*L 
, FOP li=0 TO 4 

OK Ii )=01pct*0( Ii ) 

Odp( Ii ) = (Q( Ii >-01( Ii ) )/As 
Htube( la )=Odp( Ii )/Thetab( Ii ) 

NEXT Ii 
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5060 
5064 1 
3068 1 
3072 
3076 
3080 
3084 
3088' 
3092 ‘ 



PRINTER IS 701 

RECORD TIME OF DATA TAKING 
IF ln^0 THEN 

OUTPUT 709; “TO" 

ENTER 709;Told$ 

END IF 

OUTPUT DATA TO PRINTER 



Data Set Number * " “ ,DDD»2X , 14A*‘ : J , ToldS 
Tvl Tv2 Tv3 Tldl Tld2 Tld3 
PRINT USING "10X,8( MDD.DD,2X )'• ;T(0 ),T( 1 ),T(2 ),T( 3 ) , T( 4 > , T( 5 ) , Tvav , T lav 
jbe Wall Temperatures (Deg C) Tnave Qdp 

1 2 3 4 5 6 (Deg C) (U/m'* 



3096 


PRINTER IS 701 


3100 


PRINT 


3104 


PRINT USING “10X 


5108 


PRINT 


3112 


PRINT USING "10X 


Tvav 


Tldav 


3116 


PRINT USING "10X 


3120 


PRINT 


3124 


PRINT USING ”6X, 


H 


The tab " " *' 


3128 


PRINT USING "6X, 


2) (W/m 


"2.K) (K)“““ 


3132 


Jj=0 


3136 


FOR li=0 TO Nht- 


3140 


FOR Ji*0 TO ! 


3144 


Tp( J 1 ) = T 


3148 


Jj=Jj+l 



3152 

3156 

3160 

3164 

3168 

3172 



NEXT Ji 
Jj=I i+l 
FOR Ji*0 TO 4 

Tn( J 1 )= 1 + J 1 
NEXT Ji 

PRINT USING "6X,D,lX,7(MDD.DD),lX,2(MZ.3DE).lX,l(MDD.DD>*’;Tn(Ii ),T 



p(0 ),Tp(l ),Tp(2 ),Tp(3),Tp(4 ),Tp(5 ),Tua(Ii ),0dp( Ii ) , Htubed i ) . Thetab ( 1 1 ) 



3176 


NEXT Ii 




3180 


Ok*l 




3184 


IF Im-0 THEN 




3188 


BEEP 




3192 


INPUT “OK TO STORE THIS 


DATA ! 


3196 


END IF 




32001 


J*the counter for data sets 




3204 


IF Ok=l OR Im=l THEN J-J+1 




3208 


IF Ok=l AND lm=0 THEN OUTPUT 


@Fil< 


3212 


IF Im*l OR Ok*l THEN OUTPUT 


©Plot 


3216 


Go_on=l 




3220 


IF lm=0 THEN 




3224 


BEEP 




3228 


INPUT "WILL THERE BE ANOTHER 1 


3232 


Nrun= J 




3236 


IF Gc_on=0 THEN 3272 




3240 


IF Go_on'. '0 THEN Repeat 




3244 


ELSE 




3248 


IF J Nrun+1 THEN Repeat 




3252 


END IF 




3256 


St = l 




3260 


BEEP 




3264 


. INPUT "ARE YOU SURE YOUR READY TO 


3268 


Go_on=l 




3272 


IF St >0 THEN 3280 




3275 


IF St=0 THEN GOTO 3240 




3280 


IF lm=0 THEN 
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BEEP 

PRINT 

PRINT USING •*10X/*"NOTE: data runs were stored in file 



, 10A" ; J’l ,D2f ile$ 



329G 

5300 

3304 

530B 

3312 

331G 

3320 

3324 

3328 

3332 

333G 

3340 

3344 

334B 

3352 

i 1 e 

335G 

33G0 

3364 



ASSIGN @Filel TO * 

OUTPUT @File2;Nrun-l 
ASSIGN ©Filel TO DlfileS 

ENTER @Filel;Date$,Ldtc(»),Itt, Bop ,Nht ,Natp,Nrt,Corr 
OUTPUT @File2;Date$,Ldtc( *),Itt ,Bop,Nht,NatPtNrt , Corr 
FOR 1=1 TO Nrun-1 



ENTER @Filel ;Enf( ♦ ),Uolt( » ) , Anp( ♦ ) 



OUTPUT @File2.Emf(* ),Uolt(* ),Anp(* ) 
NEXT I 

ASSIGN @Filel TO * 

PURGE “DUMMY" 



END IF 

BEEP 

PRINT 



PRINT USING ”10X,""NOTE: X-Y pairs were stored in plot data f 



, 10A" ; J“1 ,Pf ile$ 

ASSIGN @File2 TO * 
ASSIGN @Plot TO * 
BEEP 



53GB SUBEND 
3372! 

337G I XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
3380! 

33B4 DEF FNKcu(Tcu) 

338B! OFHC COPPER 

3392 Tk=Tcu+273. 15 'C TO K 

339G! Kcu=434-.ll2*Tk ' 250-300K USE FOR R-114 @2.2 C 

3400 Kcu=433,0-.l»Tk i 200-400K USE FOR R-113 @47.5 C 

3404 RETURN Kcu 
3408 FNEND 
34121 

3416 DEF FNMu(T) 

34201 CURVE FIT OF VISCOSITY 
3424 Tk-T+273.15 »C TO K 

34281 Mu=EXP(-4.4G3G+(1011.47/Tk))M.0E-3 iR-114 170-3G0 K 

3432 Mu*. 0000134 *( 10" ( 503/( Tk-2 . 15 ) ) ) 'R113 

343G RETURN Mu 

3440 FNEND 

3444 I 

3448 DEF FNCp(T) 

3452 » CURVE FIT OF Cp 

345G Tk*T+273.I5 'C TO K 

34601 Cp = . 40188+1 .65007E-3*T^ +1 .51494E-6*TI' "2-6. 67853E-10*Tk"3 'R-114 180-400 K 

3464 Cp=( 929+1 .03*T )*. 001 iP-113 

3468 Cp=Cp*1000 

3472 RETURN Cp 

3476 FNEND 

34B01 

3484 DEF FNRho( T ) 

3488 Tk=T+273.15 iC TO K 

3492 X*l-(1 .B'Tl. /7B3.95 ) 'K TO P 

34 96 f Ri} = 36.32 + 61 . 146414*X ” ( 1 /3 >+ 16 . 418015*X+ 17 . 47683S*X ' .5+1 . 1 19828 2 
3500‘ Ro=Ro/ .062428 iR-114 

3504 Ro=l .6207479E + 3-T*(2 .2l86346 + T*2.557829lE-3 ) »R-1 13 

350B RETURN Ro 
3512 FNEND 
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3516' 

3520 DEF FNPr(T) 'GOOD FOR R-114/R-113 
3524 Pr-FNCp( T )*FNMu( T )/FNK( T ) 

3528 RETURN Pr 
5532 FNEND 
3556' 

3540 DEF FNK(T ) 

5544' TOB0 K WITH T IN C 
3548 K=.071-.000261*T 

3552 RETURN K 
3556 FNEND 
35601 

3564 DEF FNTanh(Fe) 

3568 P=EXP(Fe> 

3572 Q*EXP(-Fe) 

3576 Tanh=(P-Q )/(P+Q ) 

3580 RETURN Tanh 
3584 FNEND 
5588' 

3592 DEF FNTvsv(U) 

3596 COM /Cc/ C(7) 

3600 T*C(0) 

3604 FOR 1=1 TO 7 
3608 T = T+C(I)*U'‘I 
3612 NEXT I 
3616 RETURN T 
3620 FNEND 
3624 ! 

3628 DEF FNBeta( T ) 

3632 Rop*FNRho( T+ . 1 ) 

3636 Rom=FNRho( T- . 1 ) 

3640 Beta=-2/( Rop+Rom )» ( Rop-Ron )/ . 2 

3644 RETURN Beta 

3648 FNEND 

3652 DEF FNPoly(X ) 

3656 COM /Cply/ A( 10 , 1 0 ) , C( 1 0 ) , B( 4 ) , Nop , Iprnt , Opo , I log 
3660 X1=X 
3664 Poly*B(0) 

3668 FOR 1=1 TO Nop 

3672 IF Ilog=l THEN X1=L0G(X) 

3676 Poly*Poly+B( I )*X1"I 

3680 NEXT I 

3684 IF Ilog=l THEN Poly=EXP( Pol y ) 

3688 RETURN Poly 
3692 FNEND 
3696' 

3700'XXXXXXXXXXXXXXXXXXXXXXXXXXyXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
5704 I 

3708 5UB Foiy( Df 1 le$( * ) ,Np, Itn ) 

3712 DIM R(10 ),S(10 ),Sy(12 ),5/.(12 ) 100 ), Yy( 100 ),Xy( 17 ) 

3716 COM ''Cply/ A( 10 , 10 ) , C( 10 ) , B( 4 ) , N , Iprnt , Opo , I log 

3720 COM /Xxyy/ Xpf5),Yp(5) 

3724 FOR 1=0 TO 4 

3728 B( I )=0 

3732 NEXT I 

3736 Iin=l 
3740 BEEP 

3744 input "ENTER DATA FILE NAME" ,Df i le$( 0 ) 

3743 BEEP 

3752 INPUT "ENTER NUMBER OF X-Y PAIRS", Np 
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3756 

5760 

3764 

3758 

3772 

3776 

3780 

5784 

3788 

3792 

3796 

3800 

3804 

3808 

3812 

3816 

5820 

3824 

3828 

3832 

3836 

3840 

3844 

3848 

3852 

5856 

3860 

3864 

5868 

3872 

3876 

3880 

3884 

3888 

3892 

3896 

3900 

3904 

3908 

3912 

3916 

3920 

3924 

3928 

3932 

3936 

3940 

3944 

3948 

3952 

3956 

3960 

3964 

3968 

3972 

3976 

3980 

3984 

5988 

3992 



BEEP 

INPUT "LIKE TO EXCLUDE DATA PAIRS (1 =Y , 0=N( DEFAULT .O'? ".led 
IF Ied=l THEN 
BEEP 

INPUT "ENTER NUMBER OF PAIRS TO 8E EXCLUDED'Mpex 
END IF 

ASSIGN ©File TO Dfile$(0) 

N = 2 
BEEP 

INPUT "ENTER THE ORDER OF POLYNOMIAL (DEFAULT=2) ",N 
FOR 1-0 TO N 
Sy( I )-0 
S.-.( I )=0 
NEXT I 

IF Ied=l AND Im*l THEN 
FOR I-l TO Ipex 

ENTER @File.Xy( # > 

NEXT I 
END IF 

FOR I-l TO Np-Ipex 

ENTER @FileiXy( ♦ > 

IF Opo=0 THEN 
Y*Xy( Itn-1 ) 

X»Xy( 11+Itn) 

END IF 

IF Opo=l THEN 
Y=Xy(5+Itn ) 

X-Xy( 11+Itn ) 

END IF 

IF 0po=2 THEN 
Y-Xy(5+Hn ) 

X-Xy(Itn-l ) 

END IF 

IF Ilog-1 THEN 
X-LOG(X > 

Y»LOG( Y ) 

END IF 
Xx( I >»X 
Yy( I )-Y 
R( 0 >»Y 

Sy(0) = Sy(0) + Y 
S( 1 )=X 

Sx( 1 ) = Sx( 1 ) + X 
FOR J=1 TO N 

R( J )=R( J-1 >*X 
Sy( J )-Sy( J )+R( J ) 

NEXT J 

FOR J-2 TO N*2 

S< J ) = S( J-1 >*X 
S ( J ) - S X ( J ) + S ( J ) 

NEXT J 
NEXT I 
S X ( 0 ) = Np 
FOR 1-0 TO N 
C( I ) = Sy( I ) 

FOR J-0 TO N 

A( I, J ) = S.x(I + J ) 

NEXT J 
NEXT I 

FOR 1=0 TO N-1 
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5996 CALL Dwide( I > 

4000 CALL Subtract ( I + l ) 

4004 NEXT I 

4008 B(N)-C(N)/A(N,N) 

4012 FOR 1=0 TO N-1 

4016 B(N-1-I )=C(N-1*I ) 

4020 FOR J*0 TO I 

4024 B(N-1-I )=B(N-1-I >-A( N-l-I , N- J )*B(N-J ) 

4028 NEXT J 

4052 B(N-1-I )=B(N-1-I )/A(N-l-I,N-l-I ) 

4036 NEXT I 

4040 IF lprnt=0 THEN 

4044 PRINT USING " 1 2X /•" EXPONENT COEFFICIENT * 

4048 FOR 1=0 TO N 

4052 PRINT USING 15X . DO . 5X , MO . 7DE” ; I , B(I ) 

4056 NEXT I 

4060 PRINT " " 

4064 PRINT USING 1 2X . " “DATA POINT X Y Y( CALCULATED ) 01 

SCREPANCY 

4068 FOR 1=1 TO Np 

4072 Yc=B(0) 

4076 FOR J=1 TO N 

4080 Yc = Yc+B( J )*Xx( I )"J 

4084 NEXT J 

4088 D=Yy(I)-Yc 

4092 PRINT USING '* 15X , 3D , 4X , 4( MD . 5DE , IX ) “ ; I , Xx( I ) , Yy ( I ) , Yc , D 

4086 NEXT I 

4100 END IF 

4104 ASSIGN ©File TO * 

4108 SUBEND 
41121 

4116 SUB Divide(M) 

4120 COM /Cply/ A(10, 10 > .C( 10 ) ,B( 4 ) ,N. Iprnt ,0po, I log 
4124 FOR I=M TO N 
4128 Ao=A(I,M> 

4132 FOR J«M TO N 

4136 A( I , J )=A( I , J )/Ao 

4140 NEXT J 

4144 C(I)-C(I)/Ao 

4148 NEXT I 

4152 SUBEND 

4156! 

4160 SUB Subtract(K ) 

4164 COM /Cply/ A(1 0 , 10 ) , C ( 10 ) , B( 4 ) , N , Iprn t , Opo , I log 

4168 FOR I=K TO N 

4172 FOR J=K-1 TO N 

4176 A(I,J)=A(K-1,J)-A(I,J) 

4180 NEXT J 

4184 C( I )=C(K-1 )-C( I ) 

4188 NEXT I 
4192 SUBENO 
4196! 

4200 SUB Plxn 

4204 COM /Cply/ A( 1 0 , 1 0 ) , C(10 ) , B( 4 ) , N . Iprnt , Opo , I log 
4208 COM /Xxyy/ Xx(5),Yy(5) 

4212 Pf^INTEP IS 705 
4216 BEEP 

4220 INPUT "SELECT (0 = h/h0% sane tube , l = h( HF )/h( sn > " , Irt 
4224 BEEP 

4228 INPUT "WHICH Tsat ( 1 =6 . 7 , 0 = -2 . 2 ) ‘M sa t 



162 



4232 xmn=0 
4256 Xna.'. = 10 
4240 Xstep”2 
4244 IF lrt=0 THEN 
4248 Yi>un=0 

4252 Ynax=1.4 

4256 Ystep=.2 

4260 ELSE 

4264 Ymin=0 

4268 Ynax=15 

4272 Ystep=5 

4276 END IF 
4280 BEEP 

4284 PRINT " IN ^ SPl : IP 2500 , 2200 , 8300 , 6800 ; " 

4288 PRINT "SC 0,100,0,100;TL 2,0;" 

4292 Sfx = 100/( Xnax-Xmn ) 

4296 S f y=100/ ( Ypiax-Ym n ) 

4500 PRINT "PU 0,0 PD" 

4304 FOR Xa=Xnin TO Xnax STEP Xstep 

4308 X = ( Xa-Xmi n )»S f .X 

4312 PRINT "PA";X,",0i XT;" 

4316 NEXT Xa 

4320 PRINT "PA 100,0;PU; " 

4324 PRINT "PU PA 0,0 PD" 

4328 FOR Ya=Ynin TO Ymax STEP Ystep 
4332 Y=(Ya-Ynin)»Sfy 

4336 PRINT "PA 0, " ; Y, "YT" 

4340 NEXT Ya 

4344 PRINT "PA 0,100 TL 0 2“ 

4348 FOR Xa=Xnin TO Xmax STEP Xstep 

4352 X = (Xa-Xmin )»Sfx 

4356 PRINT "PA" ; X, " , 100; XT" 

4360 NEXT Xa 

4364 PRINT "PA 100,100 PU PA 100,0 PO" 

4368 FOR Ya=Ypun TO Ynax STEP Ystep 

4372 Y = ( Ya-Ymln )*Sfy 

4376 PRINT "PD PA 1 00 . “ , Y , " YT " 

4380 NEXT Ya 

4384 PRINT "PA 100,100 PU" 

4388 PRINT "PA 0,-2 SR 1.5,2" 

4392 FOR Xa=Xnin TO Xnax STEP Xstep 
4396 X = ( Xa-Xnin )»Sfx 

4400 PRINT "PA";X,",0; " 

4404 PRINT "CP -2,-1 :LB" ;Xa; "" 

4408 NEXT Xa 

4412 PRINT "PU PA 0,0" 

4416 FOR Ya=Ynin TO Ynax STEP Ystep 
4420 IF ABS<Ya)a.E-5 THEN Ya=0 

4424 Y = ( Ya-Ynin )»Sfy 

4428 PRINT "PA 0,";Y,"“ 

4432 PRINT "CP -4,-.25.LB";Ya. "" 

4436 NEXT Ya 

4440 Xlabel$="Oil Percent" 

4444 IF lrt=0 THEN 
4448 Ylabel$="h/h0%" 

4452 . ELSE 

4456 Ylabe 1S= “ h/hsnoot h " 

4460 ENO IF 

4464 PRINT "SR 1.5,2;PU PA 50,-10 CP " ; -LEN( X labe 1 S )/ 2 . " 0 ; LB " ; X 1 abe 1 $ 
4468 PRINT "PA -11,50 CP 0 , " : -LEN( YlabelS >/2*5/6 ; ”DI 0 , 1 . LB " ■ Y1 abe 1 $ 
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4472 

447G 

4480 

4484 

4488 

4492 

4496 

4500 

4504 

4508 

4512 

4516 

4520 

4524 

4528 

4532 

4536 

4540 

4544 

4548 

4552 

4556 

4560 

4564 

4568 

4572 

4576 

4580 

4584 

4588 

4592 

4596 

4600 

4604 

4608 

4612 

4616 

4620 

4624 

4628 

4632 

4636 

4640 

4644 

4648 

4652 

4656 

4660 

4664 

4668 

4672 

4676 

4680 

4684 

4688 

4692 

4696 

4700 

4704 

4708 



PRINT “CP 0,0'* 

lpn=0 

BEEP 

INPUT "WANT TO PLOT DATA FROM A FILE ( 1 -Y , 0=N )7 " , Okp 



lcn*0 



IF Okp*l THEN 
BEEP 

INPUT "ENTER THE NAME OF THE DATA FILE” ,D_f i ie$ 
BEEP 

INPUT "SELECT (0=LINEAR, 1 =LOG( X , Y > " , I log 
ASSIGN @File TO D_file$ 

BEEP 

INPUT "ENTER THE BEGINNING RUN 
BEEP 
INPUT 
BEEP 
INPUT 
BEEP 

PRINTER IS 1 
PRINT USING 



NUMBER" ,Md 

"ENTER THE NUMBER OF X-Y PAIRS STORED" , Npaira 
"ENTER DESIRED HEAT FLUX",Q 



"4X,' 
"4X, ■ 
"4X,' 
"4X,‘ 
"4X,‘ 
“4X/ 



PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
INPUT Syn 
PRINTER IS 705 
PRINT "PU DI" 

IF Syn*l THEN PRINT 
IF Sym*2 THEN PRINT 
IF Syn=3 THEN PRINT 
Nn=4 

IF Ilog-1 THEN Nn»l 
IF Md>l THEN 

FOR I«1 TO (Md-1 > 



Select a symbol = 

1 Star 2 Plus sign""" 
Circle 4 Square""" 
Ronbus 

Right-side-up triangle" 
Up-s 1 de-doun triangle"" 



SM»" 

'SM+" 

SMo" 



ENTER @File;Xa,Ya 
NEXT I 
END IF 
01-0 

IF I log-1 THEN O-LOG(O) 

FOR I-l TO Npairs 

ENTER @File;Xa,B(^ > 
Ya-B(0> 

FOR K=1 TO Nn 

Ya-Ya+B(K )»0''K 
NEXT K 

IF Ilog-1 THEN Ya=EXP( Ya ) 
IF llog=0 THEN Ya=Ql/Ya 
IF Irt-0 THEN 
IF Xa=0 THEN 
Yo-Ya 



Ya-1 

ELSE 

Ya-Ya/Yo 
END IF 
ELSE 

H 5 m=FNHsnoot h( Q , Xa , I sa t ) 
Ya-Ya/Hsm 
END IF 



X.-. (I-l ) = Xa 
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SM" 

SR 1 . 4 , 2 . 4 ’■ 



47i: 

471B 

4720 

4724 

4729 

4732 

473S 

4740 

4744 

4749 

4752 

4756 

4760 

4764 

4768 

4772 

4776 

4780 

4784 

4788 

4792 

4796 

4800 

4804 

4808 

4812 

4816 

4820 

4824 

4828 

4832 

4836 

4840 

4844 

4848 

4852 

4856 

4860 

4864 

4868 

4872 

4876 

4880 

4884 

4888 

4892 

4896 

4900 

4904 

4908 

4912 

4916 

4920 

4924 

4929 

4932 

4936 

4940 

4944 

4948 



Vy( I-l ) = Ya 
X = ( Xa-Xfvun )#Sf/ 



Y*( Ya-Ynin )»Sfy 
IF Syn>5 THEN PRINT 
IF Syn^4 THEN PRINT 
PRINT "PA",X,Y,"" 

IF Syn>3 THEN PRINT “SR 1 . 2 , 1 . 6 " 

IF Syn=4 THEN PRINT " UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 , 8 , 4 , 0 , . " 



IF Syn*5 THEN PRINT " UC3 , 0 , 99 , -3 , -6 . -3 . 6 . 3 , 6 , 3 . -6 ; 
IF Syn*B THEN PRINT “ UC0 , 5 . 3 , 99 , 3 , ~8 , -6 , 0 , 3 , 8 . " 

IF Syn*7 THEN PRINT “UC0 , -5 . 3 , 99 . -3 , 8 , 6 , 0 , -5 , «8 j ** 

NEXT I 

BEEP 

ASSIGN @File TO * 

END IF 

PRINT “PU SM“ 

BEEP 

INPUT “WANT TO PLOT A POLYNOMIAL ( 1 = Y , 0=N , Okp 

IF Okp-1 THEN 
BEEP 

INPUT “SELECT ( 0=LINEAR . l»LOG< X , Y ) ) 'M log 
Iprnt =*1 



CALL Poly(Itn) 

FOR Xa^Xmin TO Xnax STEP Xstep/25 
Icn* Icn+ 1 



Ya'FNPoly(Xa) 
Y=( Ya-Ymin )»S f y 
X = ( Xa“Xnin )»Sf x 



IF Y<0 THEN Y*0 
IF Y>100 THEN GOTO 4868 



Pu*0 

IF Ipn=l THEN Idf=Icn MOD 2 
IF Ipn»2 THEN Idf*Icn MOO 4 
IF Ipn*3 THEN Idf*Icn MOD 8 

IF Ipn«4 THEN Idf=Icn MOD 16 

IF Ipn«5 THEN Idf^Icn MOD 32 

IF Idf-1 THEN Pu-1 

IF Pu=0 THEN PRINT “PA“ , X , Y , “PO“ 
IF Pu»l THEN PRINT “PA" , X , Y , “PU" 
NEXT Xa 
PRINT "PU" 

Ipn=Ipn+l 
GOTO 4480 
END IF 
BEEP 

INPUT “WANT TO QUIT ( 1*Y , 0-N )7 qu 1 1 
IF Iquit=l THEN 4904 
GOTO 4480 
PRINT "PU SP0“ 

SUBEND 
SUB Stats 
PRINTER IS 701 



J=0 



K=0 

BEEP 

INPUT "PLOT FILE TO ANALYZE*? ", F i ie$ 

ASSIGN @File TO FiieS 

BEEP 

INPUT “LAST RUN No*? ( 0=OU I T ) " , Nn 
IF Nn=0 THEN 5092 
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4952 Nn*Nn-J 

4956 5x=0 

49GC Sy=0 

4964 5z=0 

4968 Sx5=0 

4972 Sys=0 

4976 5zs-0 

4980 FOR 1=1 TO Nn 

4984 J=J+1 

4988 ENTER @File;Q.T 

4992 H=Q/T 

4996 5x=Sx+0 

5000 Sx5=Sxs+Q*2 

5004 Sy=Sy+T 

5008 Sys = Sys + T''2 

5012 Sz=Sz+H 

5016 Sz5=Sz5+H"2 

5020 NEXT I 

5024 Oave=5x/Nn 

5028 Tave=Sy/Nn 

5032 Have=5z/Nn 

5036 5devq = SQR( ABS( ( Nn*5xs-Sx''2 )/(Nn*(Nn-l ) ) > ) 

5040 Sdevt=SQR( ABS( (Nn*Sys-Sy"2 )/ ( Nrr< Nn-1 ) ) ) ) 

5044 Sdevh=S0R( ABS( (Nn*5zs-Sz"2 )/( Nn*( Nn- 1) ) ) ) 

5048 Sh=100*Sdevh/Have 
5052 Sq=100*Sdevq/Qave 
5056 St=100*Sdevt /Tave 
5060 IF K=1 THEN 5084 
5064 PRINT 

5068 PRINT USING •■11X,""DATA FILE: ” M4A" ; Fi le$ 

5072 PRINT 

5076 PRINT USING ■'11X,“''RUN Htube SdevH Qdp SdevQ 

5080 K=1 

5084 PRINT USING “ 1 IX ,DD , 2( 2X , D . 3DE , IX , 3D . 2D ) , 2X , DO . 3D , 1 X , 3D . 2D 

, Sq ,Tave , St 

5088 GOTO 4940 

5092 ASSIGN @Filel TO * 

5096 PRINTER IS 1 
5100 SUBEND 
5104 SUB Coef 

5108 COM /Cply/ A( 10 , 10 ) . C( 10 ) , B( 4 ) , N . Iprnt , Opo . I log 
5112 BEEP 

5116 INPUT "GIUE A NAME FOR CROSS-PLOT FILE",Cpf$ 

5120 BEEP 

5124 INPUT "OUTPUT TYPE ( 0=q vs Dt , l=h vs Dt, 2=h vs q)",Opo 
5128 CREATE BOAT Cpf$,6 
5132 ASSIGN ©File TO Cpf$ 

5136 BEEP 

5140 INPUT "SELECT ( 0 = L INEAR . 1=L0G( X , Y ) ) " , I log 
5144 BEEP 

5148 INPUT "ENTER OIL PERCENT ( - l = STOP ) " , Bop 
5152 BEEP 

5156 INPUT "ENTER TUBE NUMBER (1, 2, 3. 4, OR 5>*',Itn 
5160 IF BopK0 THEN 5176 
5164 CALL Poly( Itn ) 

5168 OUTPUT @File.Bop,B( * ) 

5172 GOTO 5144 

5176 ASSIGN @File TO ♦ 

5180 SUBEND 



Thetab SdevT"" 
J , Have , Sh , Gave 
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5184 I 



5188'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



51921 

5196 

5200 

5204 

5208 

5212 

5216 

5220 

5224 

5228 

5232 

5236 

5240 

5244 

5248 

5252 

5256 

5260 

5264 

5268 

5272 

5276 

5280 

5284 

5288 

5292 

5296 

5300 

5304 

5308 

5312 

5316 

5320 

5324 

5328 

5332 

5336 

5340 

5344 



SUB Plot 

COM /Cply/ A( 10, 10 ) ,C( 10 ) ,B(4 ),Nop, Iprnt ,0po» Hog 
DIM Xy( 17 > 

INTEGER Ii 
PRINTER IS 1 
BEEP 
I dv= 1 

INPUT "LIKE DEFAULT VALUES FOR PLOT ( 1 =Y ( DEFAULT ) , 0 = N ) 7 " . I dv 

Opo=0 

BEEP 

PRINT USING "4X, ""Select Option:""" 

PRINT USING "6X,""0 q versus de 1 1 a-T ( DEFAULT ) 

PRINT USING "6X,"”1 h versus delta-T 

PRINT USING "6X,""2 h versus q * 

INPUT Opo 
BEEP 

INPUT "SELECT UNITS ( 0*S I ( DEFAULT ), 1 -ENGL ISH )", I un 
PRINTER IS 705 
IF IdvOl THEN 
BEEP 

INPUT "ENTER NUMBER OF CYCLES FOR X-AXIS", C^ 

BEEP 

INPUT "ENTER NUMBER OF CYCLES FOR Y-AXIS", Cy 
BEEP 

INPUT "ENTER MIN X-UALUE (MULTIPLE OF 10)",Xnin 
BEEP 

INPUT "ENTER MIN Y-UALUE (MULTIPLE OF 10)".Ynin 
ELSE 

IF Opo-0 THEN 
Cy*2 
Cx-2 
Xni n» 1 
Ynin* 1000 
END IF 

IF Opo*l THEN 
Cy=2 
Cx*2 
Xn i n* 1 



5348 Ymin=100 

5352 END IF 

5356 IF Opo=2 THEN 

5360 Cy*2 

5364 C:<*2 

5368 Xrnin*1000 

5372 Ynin=100 

5376 END IF 

5330 END IF 
5384 BEEP 

5388 PRINT " IN : SP 1 ; IP 2300 , 2200 , 8300 , 6800 ; " 

5392 PRINT "SC 0, 100, 0 , 100 ; TL 2,0." 

5396 Sfx = 100/C.‘. 

5400 Sfy=100/Cy 
5404 BEEP 

5408 INPUT "WANT TO BY-PASS CAGE (1=Y, 0=NO( DEFAULT )", Ibp 
5412 IF Ibp=l THEN 5908 
5416 PRINT "PU 0,0 PD“ 

5420 Nn=9 
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5424 

5428 

5432 

545G 

5440 

5444 

5448 

5452 

5456 

5460 

5464 

5468 

5472 

5476 

5480 

5484 

5488 

5492 

5496 

5500 

5504 

5508 

5512 

5516 

5520 

5524 

5528 

5532 

5536 

5540 

5544 

5548 

5552 

5556 

5560 

5564 

5568 

5572 

5576 

5580 

5584 

5588 

5592 

5596 

5600 

5604 

5608 

5612 

5616 

5620 

5624 

5628 

5632 

5656 

5640 

5644 

5648 

5652 

5656 

5660 



FOR 1=1 TO C-+1 

Xat=Xmn*10"( I-l ) 

IF I=Cx+l THEN Nn-1 
FOR J=1 TO Nn 

IF J=1 THEN PRINT “ TL 2 0" 

IF J=2 THEN PRINT "TL 1 0" 

Xa=Xat * J 

X = LGT(Xa/Xnin 

PRINT "PA" ;X, " ,0; XT: " 

NEXT J 
NEXT I 

PRINT "PA 100,0:PU; " 

PRINT "PU PA 0,0 PD" 

Nn = 9 

FOR 1=1 TO Cy+1 

Yat=Ymn*10''< I-l > 

IF I=Cy+l THEN Nn=l 
FOR J=1 TO Nn 

IF J*1 THEN PRINT " TL 2 0" 

IF J*2 THEN PRINT "TL 1 0" 

Ya*Yat»J 

Y=LGT (Ya/Ymin )*5fy 
PRINT "PA 0,";Y,"YT" 

NEXT J 
NEXT I 

PRINT "PA 0,100 TL 0 2" 

Nn = 9 

FOR I«1 TO Cx+1 

Xat=X(nin»10''(I-l ) 

IF I=Cx+l THEN Nn=l 
FOR J*1 TO Nn 

IF J=l THEN PRINT "TL 0 2" 

IF J>1 THEN PRINT "TL 0 1“ 

Xa*Xat ♦ J 

X=LGT(Xa/Xnin)«5fx 
PRINT "PA";X,",100i XT“ 

NEXT J 
NEXT I 

PRINT "PA 100,100 PU PA 100,0 PD“ 

Nn=9 

FOR 1=1 TO Cy+1 

Yat=Ymin»10"< I-l ) 

IF I=Cy+l THEN Nn*l 
FOR J=1 TO Nn 

IF J*1 THEN PRINT "TL 0 2" 

IF J>1 THEN PRINT " TL 0 1" 

Ya=Yat*J 

Y=LGT(Ya/Ynin )*5fy 
PRINT "PD PA 100, " , Y, "YT" 

NEXT J 
NEXT I 

PRINT "PA 100,100 PU" 

PRINT "PA 0,-2 SR 1.5,2" 

I i=LGT ( Xnin ) 

FOR 1 = 1 TO C.. + 1 

Xa=Xn 1 n* 10 “ ( I-l ) 

X=LGT(Xa/Xnin)*SfA 
PRINT "PA":X,",0; " 

IF li.=0 THEN PRINT "CP -2 , -2 ; LB 10 ; PR - 
IF li>0 THEN PRINT "CP -2 , -2 . LBl 0 ; PR 0, 
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rj rj 



5664 Ii=Ii+l 

5668 NEXT I 

5672 PRINT "PU PA 0.0" 

5676 Ii-LGT(Ynin) 

5680 Y10=10 

5684 FOR 1=1 TO Cy+1 

5688 Ya=Ynin*10"(I-l > 

5692 Y=LGT(Ya/Ymn)*Sfy 

5696 PRINT "PA 0 , " : Y " 

5700 PRINT "CP -4,-.25;LB10;PR -2 , 2 : LB" ; 1 1 . " " 

5704 Ii=Ii+l 

5708 NEXT I 
5712 BEEP 
5716 Idl=l 

5720 INPUT "WANT USE DEFAULT LABELS ( 1=Y( DEFAULT ), 0=N )7 I dl 
5724 IF IdlOl THEN 
5728 BEEP 

5732 INPUT "ENTER X-LABEL " , X 1 abe 1$ 

5736 BEEP 

5740 INPUT "ENTER Y-LABEL " , Y 1 abe 1$ 

5744 END IF 
5748 IF 0po<2 THEN 

5752 PRINT "SR 1,2;PU PA 40,-14;" 

5756 PRINT "LB(T:PR -1.6,3 PD PR 1.2,0 PUiPR . 5 , -4 ; LBwo ; PR .5,1;" 

5760 PRINT "LB-T;PR . 5 , - 1 ; LBsa t ; PR .5,1;" 

5764 IF lun=0 THEN 

5768 PRINT "LB) / (K >" 

5772 ELSE 

5776 PRINT "LB) / (F)" 

5780 END IF 

5784 END IF 
5788 IF 0po=2 THEN 
5792 IF lun=0 THEN 

5796 PRINT "SR 1.5,2;PU PA 40,-14iLBq / (W/m;SR 1,1.5;PR 0.5,1;LB2;SR 1 

.5,2;PR 0.5,-l;LB)" 

5800 ELSE 

5804 PRINT "SR 1.5,2;PU PA 34,-14;LBq / (Btu/hr;PR .5,.5;LB.;PR .5, -.5; 

5808 PRINT "LBft;PR .5,1;SR 1,1.5;LB2;SR 1.5,2;PR .5,-l;LB);" 

5812 END IF 

5816 END IF 
5820 IF Opo=0 THEN 
5824 IF Iun-0 THEN 

5S28 PRINT -SR 1.5,2;PU PA -12,40.DI 0,l;LBq / (W/n;PR -1,0.5;SR 1,1. 5;L 

B2.SR 1.5, 2; PR 1, .5;LB)" 

5852 ELSE 

5836 PRINT "SR 1.5,2;PU PA -12,52;DI 0,l;LBq / (Btu/hr;PR - . 5 , . 5 : LB . ; PR 

.5, .5; " 

5840 PRINT "LBft;SR 1,1.5;PR -1,.5;LB2;PR 1,.5.SR 1.5,2;LB)" 

5844 END IF 

5848 END IF 
5852 IF Opo '0 THEN 
5856 IF lun=0 THEN 

5860 PRINT "SR 1.5,2;PU PA -12,58;DI 0,l:LBh / (U/n;PR -1,.5.SR 1,1.5;LB 

2; SR 1 .5, 2; PR .5, .5i " 

5864 * PRINT "LB.iPR .5,0:LBK)" 

5868 ELSE 

5872 PRINT "SR 1.5,2;PU PA -12,28;DI 0,l:LBh / «:Btu/hr;PP - . 5 , . 5 ; LB . ; PR 

. 5 , . 5 ; " 

5876 PRINT "LBft;PR -1,.5:SR 1,1,5;LB2;SR 1.5,2;PP .5,.5;LB.;PR .5, .5; 
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LBF 

5830 END IF 

5884 END IF 
5888 IF ldi=0 THEN 

5892 PRINT “SR 1.5.2.PU PA 50,-16 CP“ ; -LEN( X labe 1$ )/2 j " 0 . LB " ; X labe 1$ ; 

5896 PRINT “PA -14,50 CP 0, " ; -LEN( VlabelS )/2*5/6 ; "DI 0, 1 ; LB“ ; Yladel$i 

5900 PRINT "CP 0,0 DI" 

5904 END IF 
5908 lpn=0 
5912 Repeat = i 
5916 Xll=l.E+6 
5920 Xul=-l.E+6 



5924 

5928 

5932 

5936 

5940 

5944 

5948 

5952 

5956 

5960 

5964 

5968 

5972 

5976 

5980 

5984 

5988 

5992 

5996 

6000 

6004 

6008 

6012 

6016 

6020 

6024 

6028 

6032 

6036 

6040 

6044 

6048 

6052 

6056 

6060 

6064 

6068 

6072 

6077 

6080 

6084 

6068 

6092 

6096 

6100 

6104 

6108 

6112 



lcn*0 



BEEP 
Ok = l 

input “WANT TO PLOT DATA FROM A FILE ( 1 =Y( DEFAULT ), 0=N )7", Ok 

IF Ok*l THEN 

BEEP 

INPUT "ENTER THE NAME OF THE DATA FILE " , Df i le$( 0 > 

ASSIGN !?File TO Dfile$(0> 

BEEP 

Npairs=20 

INPUT “ENTER THE NUMBER OF X-Y PAIRS STORED< DEFAULT = 20 ) “ , Npairs 
BEEP 



3, 4, OR 5)“,Itn 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

INPUT 



“4X,' 
"6X,‘ 
“6X, • 
“6X, • 
“6X, ’ 
"6X, ■ 



"Select 
"1 
"3 



"5 

•6 

"7 



I t n=I tn+1 

INPUT “ENTER TUBE NUMBER <1, 

BEEP- 

PRINTER IS 1 

INPUT “WANT DEFAULT SYMBOLS? ( YES*0 ( DEFAULT ), NO* 1 )“, Synb 
Syn*I t n+2 
IF Synb*0 THEN 
GOTO 6036 
END IF 
PRINT USING 
USING 
USING 
USING 
USING 
USING 
Syn 

PRINTER IS 705 
PRINT ”PU DI" 

IF Sym*l THEN PRINT 
IF Syn*2 THEN PRINT 
IF Syn=3 THEN PRINT 
FOR 1=1 TO Npairs 
ENTER @FileiXy( * ) 

IF Opo=0 THEN 
Ya = Xy ( I tn- 1 ) 

Xa = Vy( 11 + Itn ) 

END IF 

IF Opo=l THEN 
Ya = Xy( 5+ 1 1 n ) 

Xa = Xy( 11 + Itn ) 

END IF 

IF‘Opo = 2 THEN 
Ya=Xy(5+Itn) 

Xa=Xy( I tn-1 ) 

END IF 

IF Xa^Xll THEN Xll=Xa 



symbol • 

Star 2 Plus sign""" 

Circle 4 Square 

Ronbus " " " 

Right-side-up triangle’ 

Up-s 1 de-doun triangle 



"SM*" 

"SM+" 

"SMo” 
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Bllb 

6120 

6124 

6128 

6132 

6136 

6140 

6144 

6148 

6152 

6156 

6160 

6164 

6168 

6172 

6176 

6180 

6184 

6188 

6192 

6196 

6200 

6204 

6208 

6212 

6216 

6220 

6224 

6228 

6232 

6236 

6240 

6244 

6248 

6252 

6256 

6260 

6264 

6268 

6272 

6276 

6280 

6284 

6288 

6292 

6296 

6300 

6304 

6308 

6312 

6316 

6520 

6324 

6328 

6332 

6336 

6340 

6344 

6348 

6352 



IF Xa Xul THEN Xal=Xa 
IF Iun=l THEN 

IF Opo<2 THEN Xa=Xa*l .8 
IF Opo>0 THEN Ya=Ya*.17Bl 
IF Opo*0 THEN Ya=Ya*.3l7 
IF Opo»2 THEN Xa=Xa*.317 
END IF 

X«LGT (Xa/Xpun)*Sfx 
Y-LGT( Ya/Ynin )*Sfy 
Kj=0 

CALL Synb(X,Y.Syn,Icl,Kj ) 

GOTO 6212 

IF Syn>3 THEN PRINT •’SM“ 

IF Syn<4 THEN PRINT "SR 1. 4,2.4" 

IF Icl-0 THEN 

PRINT "PA",X,Y, "" 

ELSE 

PRINT "PA“ ,X,Y, "PD" 

END IF 

IF Syn>3 THEN PRINT "SR 1.2, 1.6" 

IF Syn*4 THEN PRINT " UC2 , 4 , 99 , 0 , -8 , -4 , 0 , 0 . 8 , 4 , 0 ; " 

IF Syn*5 THEN PRINT "UC3 , 0 , 99 , -3 , -6 , -3 , 6 , 3 , 6 , 3 , -6 ; " 
IF Syn=6 THEN PRINT " UC0 , 5 . 3 , 99 , 3 , -8 , -6 , 0 , 3 , 8 i " 

IF Syn=7 THEN PRINT " UC0 , -5 . 3 , 99 , -3 , 8 , 6 , 0 , -3 , -8 ; " 
NEXT I 
PRINT "PU“ 

BEEP 



Ilab-1 

INPUT "WANT TO LABEL? ( 1*Y( DEFAULT ) ,0*N )“, I lab 
IF Ilab«l THEN 
PRINT "SP0;SP2" 

BEEP 

IF Klab«0 THEN 
Xlab=65 
Ylab=85 

INPUT "ENTER INITIAL X,Y LOCATIONS ", X lab , Ylab 

Xtt=Xlab-5 

Ytt»Ylab+8 

PRINT "SR 1,1.5- 

PRINT "SM;PA",Xtt ,Ytt,"LB Tube X File 
Ytt=Ytt-3 

PRINT "PA" ,Xtt , Ytt , "LB No Oil Name" 

IF Sym*l THEN PRINT "SM*" 

IF Sym=2 THEN PRINT "SM+" 

IF Sym=3 THEN PRINT "SMo" 

K lab = 1 
END IF 
Kj-1 

CALL Symb( X lab , Y lab , Sym , Ic 1 , K j ) 

PRINT "SR 1,1.5;SM" 

IF Sym^4 THEN PRINT "PR 2,0- 
PRINT "PR 2,0;LB"; Itn; "" 

BEEP 

INPUT "ENTER BOP( 0-DEFAULT >", Bop 
IF Bop<10 THEN PRINT "PR 3 , 0 ; LB " : Bop ; " " 

. IF Bop^9 THEN PRINT "PR 1 . 5 , 0 i LB" . Bop ; " " 

PRINT "PR 2,0.LB";Dfile$(0): "" 

PRINT "SP0;SP1;SR 1.5,2" 

Ylab=Ylab-5 
END IF 
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6360 

6364 

6368 

6372 

6376 

6380 

6384 

6388 

6392 

6396 

6400 

6404 

6408 

6412 

6416 

6420 

6424 

6428 

6432 

6436 

6440 

6444 

6448 

6452 

6456 

6460 

6464 

6468 

6472 

6476 

6480 

6484 

6488 

6492 

6496 

6500 

6504 

6508 

6512 

6516 

6520 

6524 

6528 

6532 

6536 

6540 

6544 

6548 

6552 

6556 

5560 

6564 

5568 

5572 

6576 

6580 

6584 

6538 

6592 



ASSIGN @File TO * 

Xll-Xll/1 .2 
Xul*=Xui*l .2 
I GOTO 8040 
END IF 

PRINT ”PU sri" 

BEEP 
Go_on= 1 

INPUT "WANT TO PLOT A POLYNOMIAL (1=Y ( DEFAULT ), 0=N )7 \ Go_on 
IF 6o_on=l THEN 
BEEP 

PRINTER IS 1 

INPUT "WANT DEFAULT LINE TYPE? (YES-0 ( DEFAULT ). NO- 1 )", Ln 
Ipn=I t n 
IF Ln-0 THEN 
GOTO 6448 
END IF 

PRINT USING "4X/’ "Select line type^ 

PRINT USING "6X,""0 Solid line""" 

PRINT USING "6X,""1 Dashed 

PRINT USING "6X,""2»,,5 Longer line - dash 

INPUT Ipn 
PRINTER IS 705 
BEEP 
Ilog=l 

INPUT "SELECT ( 0-LIN . l“LOG( DEFAULT ))", I log 
Iprnt=l 

CALL Poiy( Df 1 le$( * ) , Npa irs , I tn ) 

FOR Xx-0 TO Cx STEP Cx/200 
Xa = Xnin* 10''X^. 

IF Xa<Xli OR Xa>Xul THEN 6572 

Icn=Icn+l 

Pu=0 

IF Ipn=l THEN Idf-Icn MOD 2 
IF Ipn-2 THEN Idf-Icn MOD 4 
IF Ipn=3 THEN Idf-Icn MOD 8 

IF Ipn=4 THEN Idf-Icn MOD 16 

IF Ipn-5 THEN Idf-Icn MOD 28 

IF Idf»l THEN Pu=l 

Ya-FNPoly(Xa ) 

IF Ya^^Ynin THEN 6572 
IF Iun-1 THEN 

IF Opo<2 THEN Xa-Xa*l .8 
IF Opo>0 THEN Ya=Ya*.1761 
IF Opo=0 THEN Ya=Ya*.3l7 
IF 0po*2 THEN Xa=Xa*.317 
END IF 

Y-LGT (Va/yrnin)*Sfy 

X = LG T X a / X m 1 n ) * S f A 

IF X'^0 THEN Y=0 

IF YU00 THEN GOTO 6572 

IF P-j-0 THEN PRINT *' PA X , Y , " PD " 

IF Pu=l THEN PRINT " PA " , X , Y , " PU " 

NEXT X,. 

. PRINT "PU" 

END IF 
BEEP 

INPUT "WANT TO QUIT ( 1-Y , 0-N( DEFAULT )) ‘Mqt 
IF Iqt=l THEN 6600 
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6556 60 TO -516 

6600 PRINT "FU PA 0,0 SP0" 

6604 SUBEND 
6608' 

6612 ‘ xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 



6616* 






6620 


SUB Sy(nb(X,Y,Syn,Icl,Kj ) 


6624 


IF Syn'3 


THEN PRINT "SM" 


6628 


IF Syn<'4 


THEN PRINT "SR 1. 4,2.4" 


6632 


V ad*0 




6656 


IF Kj*l 


THEN Yad=.8 


6640 


IF lcl=0 


THEN 


6644 


PRINT 


"PA" ,X , Y+Yad, " " 


6648 


ELSE 




6652 


PRINT 


"PA" ,X, Y+Yad, "PD" 


6656 


END IF 




6660 


IF Syn/5 


THEN PRINT "SR 1.2, 1.6" 


6664 


IF Syn=4 


THEN PRINT " UC2 , 4 , 99 , 0, -8 , -4 , 0 , 0 , 8 , 4 , 0 ; " 


6668 


IF Syn=5 


THEN PRINT "UC3 , 0 , 99, -3 , -B , -3 , 6 , 3 , G , 3 . -G , 


6672 


IF Syn=6 


THEN PRINT "UC0, 5 . 3 , 99 , 3 , -8 , -G , 0, 3 , 8 ; " 


6676 


IF Syn*7 


THEN PRINT "UC0 , -5 . 3 , 99 . -3 , 8 , 6 , 0 , -3 , -8 ; " 


6680 


IF Kj*l 


THEN PRINT "SM;PR 0,-.8" 


6684 


SUBEND 




6688! 







6692 ‘XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
6696! 

6700 SUB Fixup 
6704! FILE: FIXUP 
6708 J 

6712 DIM Enf(34),Anp( 11 ),Uolt<4 ),Ldtc<4) 

6716 BEEP 

6720 INPUT "OLD FILE TO FIXUP ", D2f i le$ 

6724 ASSIGN @File2 TO D2file$ 

6728 Dlfile$="TEST" 

6732 CREATE BOAT Dlfile$,60 
6736 ASSIGN @Filel TO DlfileS 

6740 ENTER @Fi le2 ; Nrun , Da t e$ , Ldt c ( ♦ ) , 1 1 1 , Bop , Nht , Na t p , Nrt , Corr 
6744 OUTPUT @Fi 1 el ; Nrun , Dat e$ , Ld t c ( * ) , 1 1 1 , Bop , Nht , Nat p , Nrt , Corn 
6748 FOR 1*1 TO Nrun 

6752 ENTER @F 1 1 e2 ; To 1 d$ , Enf ( * > , Uol t < * ) , Amp< ♦ ) 

6756 IF 1*1 THEN 6764 

6760 OUTPUT @Filel .Bop,Told$,Enf( * >,Uolt( * >,Amp( * ) 

6764 NEXT I 

6768 ASSIGN @File2 TO * 

6772 ASSIGN @Filel TO * 

6776 ‘ RENAME "TEST" TO D2_file$ 

6780 BEEP 2000, .2 
6784 BEEP 4000, .2 
6788 BEEP 4000, .2 
6792 SUBEND 

6796 I 

5800 i^XYyYXXVXKXxXXXXXXXXXXXXXyxyyxyXXyyxXXXXXXXKvyyyyHVXxyyYxxyyXXxXXyyv 
6304 I 

6808 SUB Move 
68121 FILE NAME MOUE 
6816' 

6320 DIM A'' 66 / ,B(66 ^,C^66 ),D(B6 ),E(66 ),F‘. 66 ),G(66 >,H<,66 ),J(66 66 ),L^66 ^,M(66 

) 

6824 DIM N(66 »,Enf^ 34 ),Uolt( 2 ),Amp( 11 ,Ldtc(4) 

6528 BEEP 
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6852 INPUT "OLD FILE TO MOUE ' , D2_f i 1 e$ 

6836 ASSIGN @Fiie2 TO D2_fiie$ 

6840 ENTER @Fi le2 ; Nrun, ToldS ,Ldtc( * > , 1 1 1 , Bop , Nht ,Natp ,Nrt ,Corr 

6844 FOR I-i TO Nrun 

6848 ENTER @File2;Told$ 

6852 ENTER @Fi i e2 ; A< I ) , B( I ) . C< I ) , D( I ) , E( I ) , F( I ) , G( I > . H( I ) , J ( I ) , K ( I ) , LU ) , M( 

I ).N< I ) 

6856 ENTER @Fi le2 ; Enf ( * ) , Uo 1 1 ( * ) , Anp( ♦ ) 

6860 NEXT I 

6864 ASSIGN @File2 TO * 

6868 BEEP 

6872 INPUT "SHIFT DISK AND HIT CONTINUE",Ok 
6876 BEEP 

6880 INPUT "INPUT BOAT SIZE'\Size 
6884 CREATE BOAT D2_f i le$ , Size 

6838 ASSIGN @Filel TO D2_fiie$ 

6892 OUTPUT @Fi lei ; Nrun ,Date$,Ldtc( * ) , I t t ,Bop , Nht , Natp , Nrt ,Corr 

6896 FOR I-l TO Nrun 

6900 OUTPUT @Filel.Told$ 

6904 OUTPUT @Filel;A( I ),B( I ),C( I ),0( I ),E( I ),F( I ),G( I ),H( I ), J( I >,K( I ),L( I ),M 

( I >,N( I ) 

6908 OUTPUT @Filel;Emf( » ),Uolt< * ).Afnp( * ) 

6912 NEXT I 

6916 ASSIGN @Filel TO * 

69201 RENAME "TEST" TO D2_fileS 

6924 BEEP 2000, .2 

6928 BEEP 4000 2 

6932 BEEP 4000,. 2 

6936 PRINT "DATA FILE MOUED" 

6940 SUBEND 
69441 

69481XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

6952! 

6956 SUB Purg 
6960 BEEP 

6964 INPUT "ENTER FILE NAME TO BE DELETED" , Fi le$ 

6968 PURGE File$ 

6972 GOTO 6960 
6976 SUBEND 

69801 

69841 XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
6988! 

6992 SUB Comb 
6996! FILE NAME: COMB 
7000 ‘ 

7004 DIM Emf(34 ),Uolt(2 ),Amp< 11 ),Ldtc(4) 

7008 BEEP 

7012 INPUT "OLD FILE TO F I XUP " , DZ_f i 1 e$ 

7016 ASSIGN aFiie2 TO D2_file$ 

7020 Ol^f ile$="TE5T" 

7024 CREATE BOAT Dl_fiie$,30 
7028 ASSIGN (9Filel TO Dl^fileS 

7032 ENTER @Fi le2 ; Nrun , Da t e$ , Ldtc ( ♦ > , I t t , Bop , Nh t , Nat p , Nrt , Corn 

7036 IF K=0 THEN OUTPUT @Fi 1 e 1 ; Nrun , Da te$ , Ld t c ( * ) , 1 1 1 , Bop , Nh t , Na t p , Nr t , Corr 

7040 FOR 1=1 TO Nrun 

7044 ENTER @Fi le2 ; Bop , Tol d$ , Emf < * ) , Uo 1 1 ( ♦ ) , Amp ( ♦ ) 

7048 OUTPUT @Filel : Bop ,ToldS,Emf ( * ),Uclt( * ),Amp( * ) 

7052 NEXT I 

7056 ASSIGN @File2 TO * 

70601 RENAME "TEST" TO D2_file$ 
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70G4 BEEP 4000, 

70G8 BEEP 
7072 Oka=l 

707B INPUT "WANT TO ADD ANOTHER FILE (1*Y , 0=N( def aul t ) )7 “ , Oka 
7080 IF Oka=l THEN 
7084 K*1 

7088 BEEP 

7092 INPUT "GIUE NEW FILE NAME",Nfile$ 

7096 ASSIGN @File2 TO NfileS 
7100 GOTO 7052 
7104 END IF 

7108 ASSIGN @File2 TO * 

7112 SUBEND 
711G‘ 

7120 f XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



7124 ! 






7128 


SUB 


Readplot 


7132 


DIM 


0dp(5 ),Htube(5),Thetab(S) 


713G 


PRINTER IS 701 


7140 


INPUT "ENTER FILE NAME",FileS 


7144 


INPUT "ENTER THE NUMBER OF DATA PAIRS”, Nrun 


7148 


ASSIGN @Filel TO File$ 


7152 


FOR 


1=1 TO Nrun 


7156 




ENTER @Filel :0dp< • ),Htube( • ),Thetab(* ) 


7160 




PRINT Odp( * ) 


7164 




PRINT 


7168 




PRINT Htube( * ) 


7172 




PRINT 


7176 




PRINT Thetab( * ) 


7180 




PRINT 


7184 




PRINT 


7188 


NEXT I 


7192 


SU8END 
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